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ABSTRACT 

In order to generate credible 0.1-2/im SEDs, the GAMA project requires many Giga- 
bytes of imaging data from a number of instruments to be re-processed into a standard 
format. In this paper we discuss the software infrastructure we use, and create self- 
consistent ugrizYJHK photometry for all sources within the GAMA sample. Using 
UKIDSS and SDSS archive data, we outline the pre-processing necessary to standard- 
ise all images to a common zeropoint, the steps taken to correct for seeing bias across 
the dataset, and the creation of Gigapixel-scale mosaics of the three 4x12 deg GAMA 
regions in each filter. From these mosaics, we extract source catalogues for the GAMA 
regions using elliptical Kron and Petrosian matched apertures. We also calculate Sersic 
magnitudes for all galaxies within the GAMA sample using SIGMA, a galaxy component 
modelling wrapper for GALFIT 3. We compare the resultant photometry directly, 
and also calculate the r band galaxy LF for all photometric datasets to highlight the 
uncertainty introduced by the photometric method. We find that (1) Changing the ob- 
ject detection threshold has a minor effect on the best-fitting Schechter parameters of 
the overall population (M* ±0.055 mag, a ±0.014, (j>* ±0.0005 h 3 Mpc~ 3 ). (2) An offset 
between datasets that use Kron or Petrosian photometry regardless of the filter. (3) 
The decision to use circular or elliptical apertures causes an offset in M* of 0.20 mag. 

(d) TVip hp^t.-fittince f"\ptipplitpr nflraxnpt.prc; from tnt.Pil-mpicmit.iirlp nVint.nmpt.rip svstpms 
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1 INTRODUCTION 

When calculating any statistic it is essential that the sample 
used to generate it is both numerous and without systematic 
bias. For a number of fundamental parameters in cosmology, 
for example the galaxy stellar mass function or the total lu- 
minosity density, the dataset used will be made up of a large 
sample of galaxi es, and contain a measure of the flux from 
each galaxy (e.g.. lHill et allioiOT ). Unfortunately, our ability 
to accurately calculate the flux of any galaxy is imprecise; at 
some distance from its centre the luminosity of the galaxy 
will drop into the background noise and the quantification 
of the missing light beyond that point is problematic with 
no obviously correct proce dure. Even using deep photome- 
try (fiB > 29 mag arcsec 2 ) . [Caon et ah! (|l990T l did not reveal 
the presence of an elliptical galaxy light profile truncation. 
A number of methods to calculate the flux from a galaxy 
have been proposed. They tend to be either simple and im- 
practical, such as setting the aperture to be a fixed constant 
size, or limit it using a detection threshold (ignoring the 
missing light issue completely), or complex and subject to 
bias, such as using the light distribution of the easily de- 
tected part of th e object to calculate the si ze the aperture 
should be set to (|Petrosianlll976l ; lKronlll980r ) . which will re- 
turn a different fraction of the total light emissio n depend- 
ing on whether the g ala xy follows an exponent ial (jPattersonl 



1940; Freeman 



the galaxy tollows an exponential (Fattcrsoi 
197(1) o r lde~Vaucouleurs) (|l948l ) light profile 



Cross fc Driver ( 20021 ') discuss the use of different missing- 
light estimators and their inherent selection effects. A third 
option is to attempt to fit a light profile, such as the afore- 
mentioned deVauc ouleur and exponential profiles (i.e. SDSS 
model magnitudes, Stoughton et alj 2002 ), or the more gen- 
eral Sersic profile (|Sersidll963l . iGraham fc Driver! l2005h . to 
the available data, and integrate that profil e to infinity to 
calcul ate a total-magnitude for the galaxy. IGraham et al.l 
(2005) investigate the discrepancy between the Sersic and 
SDSS Petrosian magnitudes for different light profiles, pro- 
viding a simple correction for SDSS data. 
Unfortunately, no standard, efficacious photometric formula 
is used in all surveys. If one looks at 3 of the largest pho- 
tometric surveys, 2MASS, SDSS and UKIDSS, one finds 
a variety of methods. The 2MASS survey dataset contains 
Isophotal and Kron circularised, elliptical aperture magni- 
tudes (elliptical apertures with a fixed minimum semi-minor 
axis), and an elliptical Sersic total magnitude. SDSS use two 
methods for their extended source photometry: petromag, 
which fits a circular Petrosian aperture to an object, and 
modelmag, which chooses whether an exponential or de- 
Vaucouleur profile is the more accurate fit and returns a 
magnitude determined by integrating the chosen profile to 
a specified number of effective radii (profiles are smoothly 
truncated between 7 and 8 R e for a deVaucouleur profile, 
between 3 to 4 R e for an exponential profile). The mod- 
ELMAGS used within this paper specify the profile type in 
the r band, and use that profile in each passband. UKIDSS 
catalogues were designed to have multiple methods: again 
a circular Petrosian magnitude and a 2D Sersic magnitude, 
calculated by fitting the best Sersic profile to the source. 
The 2D Sersic magnitude has not yet been implemented. As 
these surveys then form the basis for photometric calibra- 
tion of other studies it is important to understand any biases 
that may be introduced by the photometric method. 



The GAMA survey (|Driver et all l2009h is a multi- 
wavelength (151.6 nm— ~ 6m), spectroscopic survey of 
galaxies within three 4 x 12 deg regions of equatorial sky 
centred around 9h, 12h and 14. 5h (with aspirations to es- 
tablish further blocks in the SGP). Amongst other legacy 
goals, the survey team will create a complete, magnitude 
limited sample of galaxies with redshift and colour infor- 
mation from the FUV to Radio passbands, in order to accu- 
rately model the AGN, stellar, dust and gas contents of each 
individual galaxy. This requires the combination of obser- 
vations from many surveys, each with different instrument 
resolutions, observational conditions and detection technolo- 
gies. As the luminous output of different stellar populations 
peaks in different parts of the EM spectrum, it is not always 
a simple task to match an extended source across surveys. 
SDSS, which covers only a relatively modest wavelength 
range (300 — 900 nm), detects objects using a combination 
of all filters, defines apertures using the r band and then 
applies them to ugiz observations to negate this problem. 
This ensures a consistent deblending outcome and accurate 
colours. The UKIDSS extraction pipeline generates indepen- 
dent detection lists separately in each frame (i.e., for every 
filter) and merges these lists together for frames that cover 
the same region of sky (a frame set). Sources are then de- 
fined as dete ctions within a certa in tolerance. This process 
is detailed in lHamblv et al.l ((2008) . Unfortunately, it is sus- 
ceptible to differing deblending outcomes that may produce 
less reliable colours. As a key focus of GAMA is the pro- 
duction of optimal SEDs, it is necessary for us to internally 
standardise the photometry so that is immune to aperture 
bias from u-K. The pipeline outlined in this paper is the 
result of these efforts. 

Imaging data is taken from UKIDSS DR4/SDSS DR6 ob- 
servations. In section [2] we briefly outline the surveys that 
acquired the data we use in this work. In section [3] we 
describe how we standardised our data and formed image 
mosaics for each filter/region combination. Section [4] dis- 
cusses the photometric methods we use, and in sections [S] 
and [15] we discuss the source catalogues produced following 
source extraction on these mosaics. We define the photom- 
etry we are using as the GAMA standard in section [7] Fi- 
nally, in order to quantify the systematic bias introduced 
by the choice of photometric method, we present r band 
luminosity functions, calculated using a series of different 
photometric methods, in Section [5] Throughout we adopt 
an h — 1, £Im = 0.3, Qa = 0.7 cosmological model. All mag- 
nitudes are quoted in the AB system unless otherwise stated. 
Execution speeds provided are from a run of the pipeline on 
a 16 processor server built in 2009. As other processes were 
running simultaneously, processing speed will vary and these 
parameters should only be taken as approximate timescales. 



2 SURVEY DATA 
2.1 GAMA 



The GAMA project dDriver et all l2009h aims to study 
galaxy formation and evolution using a range of cutting-edge 
instruments (A AT, VST, VISTA, ASKAP, HERSCHEL 
WISE, GALEX and GMRT), creating a database of ~350 
thousand galaxies observed from UV to radio wavelengths. 
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Figure 1. Coverage of the equatorial region of sky that contains the GAMA regions, by SDSS DR6 (blue), UKIDSS LAS (red) and 
GALEX (violet) imaging. The three rectangular boxes contain the regions of sky surveyed by GAMA. 



The first stage of the GAMA project, GAMA I, covers 
144 deg 2 of equatorial sky, in three separate 4 deg x 12 deg 
regions centred at 9h +ld (GAMA9), 12h +0d (GAMA12) 
and 14h30m +0d (GAMA15). These areas have complete 
SDSS coverage, and we are in the process of obtaining com- 
plete UKIDSS-LAS coverage (Figure Q). Between 2008 and 
2010, the GAMA project was allocated 66 nights on the 
AAT to use the AAOmega spectrograph in order to carry 
out the GAMA I spectroscopic campaign. 
A complete description of the inpu t catalogue fo r the s pec- 
troscopic campaign can be found in lBaldrv et al.l (|2010t) . To 
summarise: the aim is to provide spectroscopy of all galaxies 
in the GAMA I regions brighter than r po t r o,SDSS = 19.4 mag, 
•Zmodd.SDSS = 18.2 mag and 7f k ron,AB = 17.5 mag, with the 
sample extended to r pct ro,SDSS < 19.8 mag in the GAMA12 
region. Where a galaxy would not be selected by its r mag- 
nitude, but would be selected using the K or z magnitude 
cut, the galaxy must also have r po t ro ,SDSS < 20.5 mag. This 
ensures that the galaxy is credible and the likelihood of 
obtaining a redshift is reasonable. In order to guarantee a 
complete sample of galaxies, including compact objects, the 
GAMA input catalogue utilises a star-galaxy selection al- 
gorithm that includes optical (r ps f-T model , g-i) an d infrared 
colour selections (J-K). The latter uses colours taken from 
sources extracted using this pipeline. 

The 2008 and 2009 observations make up a sample of 100051 
reliable redshifts, of which 82696 come from the AAOmega 
spectrograph. The tiling strategy used to allocate object s 
to AAOmega fibres is detailed in iRobotham et al l (|2010h . 
A breakdown of redshift completeness by luminosity and 
colo ur selection o f the year 2 observations is shown in Table 
5 of iBaldrv et al.l . and in Table 3 of the same paper there is 
a list of spectra used from external surveys. 



2.2 SDSS 



The Sloan Digital Sky Survey (SPSS. lYork et alj|2000h is 
the largest combined photometric and spectroscopic survey 
ever undertaken, and contains spectra of 930 000 galaxies 
spread over 8423 square degrees of sky, with imaging in five 



filters with effective wavelengths between 300 and 1000 nm 
(ugriz). SDSS dat a has been publ i cly re leased in a series 
of 7 data releases. [Abazaiian et al.l (|2004 ) state that SDSS 
imaging is 95% complete to it = 22.0 mag, g = 22.2 mag, 
r — 22.2 mag, i = 21.3 mag, and z = 20.5 mag (all depths 
measured for point sources in typical seeing using the SDSS 
photometric system, which is comparable in all bands to the 
AB system ±0.05 mag). 

Images are taken using an imaging array of 30 2048x2048 
Tektronix CCDs with a pixel size of 0.396 arcsec, but only on 
nights where the seeing is < 1.5 arcsec and there is less than 
1% uncertainty in the zeropoint. When such conditions are 
not reached, spectroscopy is attempted instead. SDSS cat- 
alogues can be accessed through the SDSS Catalog Archive 
Server (CAS), and imaging data through the Data Archive 
Server (DAS). 

Astrometry for SDSS-DR6 (|Pier et al.ll2003l ) is undertaken 
by comparing r band ob servations to the USN O CCD Astro- 
graph Catalog (UCAC, IZacharias et alJl200Ch. where it had 
coverage at time of release, or Tvcho-2 <|H0g et all l200Ch , 
when UCAC did not have coverage. For sources brighter 
than r = 20 mag, the astrometric accuracy when compar- 
ing to UCAC is 45 mas, and when comparing to Tycho-2 is 
75 mas. In both cases, there is a further 30 mas systematic 
error, and a relative error between filters (i.e., in ugiz) of 
25-35 mas. 

The GAMA input catalogue is defined using data 
from the sixth data release catalogue (SDSS DR6, 
lAdelman-McCarthv et~ai1 120081 ). The GAMA regions fall 
within the SDSS DR6 area of coverage, in SDSS stripes 9 to 
12. 



2.3 UKIDSS 

UKIDSS l|Lawrence et all2007h is a seven year near-infrared 
(NIR) survey programme that will cover several thousand 
degrees of sky. The programme utilises the Wide Field 
Camera (WFCAM) on the 3.8 m United Kingdom Infra- 
Red Telescope (UKIRT). The UKIDSS program consists 
of five separate surveys, each probing to a different depth 
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and for a different scientific purpose. One of these surveys, 
the UKIDSS Large Area Survey (LAS) will cover 4000 deg 2 
and will overlap with the SDSS stripes 9 to 16, 26 to 33 
and part of stripe 82. As the GAMA survey regions are 
within SDSS stripes 9 to 12, the LAS survey will provide 
high density near-IR photometric coverage over the entire 
GAMA area. The UKIDSS-LAS survey observes to a far 
greater depth [Klas = 18.2 mag using the Vega magni- 
tude system) than the previous Two- Micron All-Sky Survey 
(2MASS, K2MASS = 15.50 mag using the Vega magnitude 
system) . 

When complete, the LAS has target depths of K = 18.2 mag, 
H = 18.6 mag, J = 19.9 mag (after two passes; this paper 
uses only the first J pass which is complete to 19.5 mag) 
and Y = 20.3 mag (all depths measured use the Vega sys- 
tem for a point source with 5a detection within a 2 arcsec 
aperture). Currently, observations have been conducted in 
the equatorial regions, and will soon cover large swathes of 
the Northern Sky. It is designed to have a seeing FWHM of 
< 1.2 arcsec, photometric rms uncertainty of < ±0.02 mag 
and astrometric rms of < ±0.1 arcsec. Each position on the 
sky will be viewed for 40s per pass. All survey data for this 

paper is taken from the fourth data rel ease (DR4). 

The WFCAM Science Archive (WSA. lHamblv et alj|200§ ) 
is the storage facility for post-pipeline, calibrated UKIDSS 
data. It provides users with access to fits images and CASU- 
generated object catalogues for all five UKIDSS surveys. We 
do not use the CASU generated catalogues for a few rea- 
sons. Firstly, the CASU catalogues for early UKIDSS data 
releases suffer from a fault where deblended objects are sig- 
nificantly brighter t han their parent object, in some ca ses by 
several mag nitudes (|Smith et al.ll2009l ; lHill et al.ll2010l ). Sec- 
ondly, the CASU catalogues only contain circular aperture 
fluxes. Thirdly, CASU decisions (e.g., deblending and aper- 
ture sizes) are not consistent between filters. For instance, 
the aperture radius and centre used to calculate kpetro- 
MAG of a source is not necessarily the same as the aperture 
radius and centre used to calculate ypetromag or hpetro- 
MAG. We require accurate extended-aperture colours; the 
CASU catalogues do not provide this. 



3 CONSTRUCTION OF THE MOSAICS FROM 
SDSS AND UKIDSS IMAGING 

One of GAMA's priorities is the accurate measurement 
of SEDs from as broad a wavelength range as possible. 
This is non-trivial when combining data from multiple sur- 
veys. While each survey may be internally consistent with 
data collected contemporaneously, conditions between sur- 
veys can vary. In particular, seeing and zeropoint parameters 
may greatly differ between frames. When matching between 
surveys one may find an object in the centre of the frame 
in one survey is split across two frames in another survey. 
There may also be variation in the angular scale of a pixel 
between different instruments, and even when two instru- 
ments have the same pixel size, a shift of half of a pixel 
between two frames can cause significant difficulties in calcu- 
lating colours for small, low surface brightness objects. Fur- 
thermore, in order to use SExtractor in dual frame mode, 
the source-detection and the observation frame pixels must 
be calibrated to the same world coordinate system. In the 
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Figure 2. A histogram of the calculated total zeropoints for the 
fields used to create our master region mosaics. 



GAMA survey, we have attempted to circumvent these dif- 
ficulties by creating Gigapixel scale mosaics with a common 
zeropoint and consistent WCS calibration. The construction 
process is outlined within this section. 

To generate our im age mosaics, we use the Terapix SWARP 
ijBertin et alJl2002h utility. This is a mosaic generation tool, 
and how we utilise it is described in subsection I3.4I Before 
we can use this software, however, it is necessary for us to 
normalise the contributing SDSS and UKIDSS data to take 
into account differences in sky conditions and exposure times 
between observations. For every file we must identify its cur- 
rent zeropoint (see the distribution in Figure [2|, and trans- 
form it to a defined standard. This process is described in 
subsections I3.ll and I3.2I 



3.1 UKIDSS: Acquisition of data and 

renormalisation to a common zeropoint 

UKIDSS imaging is stored within the WFCAM Science 
Archive (WSA). We downloaded 862 Y, 883 J, 931H 
and 928 K band compressed UKIDSS-LAS fits files that 
contained images of sky from the GAMA regions. These 
files were decompressed using the imcopy utility. The files 
for each band are stored and treated separately. 
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Band 


AB offset (mag) 


u 


-0.04 


g 





r 





i 





z 


+0.02 


Y 


0.634 


J 


0.938 


H 


1.379 


K 


1.900 



Table 1. Conversion to AB magnitudes. The SDSS photometric 
system is roughly equivalent to the AB magnitude system, with 
only small offsets in the u and z passbands. UKIDSS photometry 
is calcula ted on the Vega mag nitude system, and our conversions 
are from iHewett et alj [120061) . Whilst we convert UKIDSS data 
using a high precision parameter, it should be noted that the con- 
versio n uncertainty is only known to ~ ±0.02 mag llCohen et al.l 
l2003h . 



A specially designed pipeline accesses each file, reads 
the MAGZPT (z7P mag ), EXP.TIME (t), airmass 
(0.5 * (AMSTART + AMEND) = secxmean) and 
EXTINCT (Ext) keywords from the fits header and creates 
a total AB magnitude zeropoint for the file using Equation[T] 



ZP 



u „ !g -2.51og(-)-Extx (secxn 



-1) + ABV X 
(1) 

where ABVx is the AB magnitude of Vega in the X band 
(Tabled]). 

To correct each frame to a standard zeropoint (30), the 
value of each pixel is multiplied by a factor, calculated 
using Equation [2] Whilst we show the distribution of frame 
zeropoints in Figure [2] in bins of 0.1 mag, we use the actual 
zeropoint of each frame to calculate the total AB magnitude 
zeropoint. This has a far smaller variation ( e.g., 0.02 mag 
in ph otometric conditions in the JHK filters; rWarren et al.l 
12007ft . 



pixelmodifier = 10 °' 4 ' ZPt ' 



-30) 



(2) 



A new file is created to store the corrected pixel table, 
and the MAGZPT fits header parameter is updated. The 
SKYLEVEL and SKYNOISE parameters are then scaled 
using the same multiplying factor. This process takes 3 sec- 
onds per file. 



3.2 SDSS: Acquisition of data and 

renormalisation to a common zeropoint 

The tsField and fpC files for the 12757 SDSS fields that cover 
the GAMA regions were downloaded from the SDSS data 
archive server (das.sdss.org) for all five passbands. Again, 
the files for each passband are stored and treated separately. 

We use a specially designed pipeline that brings in the 
aa (zeropoint), kk (extinction coefficient) and airmass key- 
words from a field's tsField file, and the EXPTIME (t) key- 
word from the same field's fpC file. Combining these using 
Equation [3] we calculate the current total AB magnitude 
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Figure 3. A histogram of the calculated seeing of the fields used 
to create our master region mosaics. 



zeropoint of the field (ZPtotai)- 

ZP total = — aa — 2.51og(l/t) — kk x airmass + sAo (3) 

where sAo is the offset between the SDSS magnitude sys- 
tem and the actual AB magnitude system (—0.04 mag for u, 
0.02 mag for z, otherwise 0). The SDSS photometric zero- 
point uncert ainty is estimated to be no larger than 0.03 mag 
m any band jlvezic et al.ll2004h . We calculate the multiplier 
required to transform every pixel in the field (again using 
Equation [2]) to a standard zeropoint (30). As every pixel 
must be modified by the same factor, we utilise the f carith 
program (part of the Ftools package), to multiply the en- 
tire image by pixelmodifier. f carith can normalise an SDSS 
image every 1.25 seconds. 



3.3 Correction of seeing bias 

As observations were taken in different conditions there 
is an intrinsic seeing bias between different input images, 
and between different filters (Figure [3j. This could cause 
inaccuracies in photometric colour measurements that use 
apertures defined in one filter to derive magnitudes in a 
second filter. To rectify this problem, it is necessary for 
us to degrade the better quality images to a lower seeing. 
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However, if we degrade all images to our lowest quality 
seeing (3.12 arcsec), we should lose the ability to resolve 
the smallest galaxies in our sample. Therefore, we elect 
to degrade our normalised images to 2 arcsec seeing. The 
fraction of images with seeing worse than 2 arcsec is 4.4%, 
2.7%, 2.5%, 2.1%, 1.7%, 0%, 0%, 1.3% and 0.9% in u, g, r, 
i, z, Y, J, H and K, respectively. Images with worse seeing 
than 2 arcsec are included in our degraded seeing mosaics. 
We do not attempt to modify their seeing. Although each 
survey uses a different method of calculating the seeing 
within their data (SDSS use a double gaussian to model 
their PSF, UKIDSS use the average FWHM of the stellar 
sources within the image), we assume that the seeing 
provided for every frame is correct. 

To achieve a final PSF FWHM of 2 arcsec (af ina i) we 
assume that the seeing within an image follows a perfect 
Gaussian distribution, Oiniuai- Theoretically, a Gaussian 
distribution can be generated from the convolution of two 
Gaussian distributions. The fgauss utility (also part of 
the ftools package) can be used to convolve an input 
image with a circular Gaussian with a definable standard 
deviation (a req ), calculated using Equation [4] 



are 1 — Y "jinal a initial (4) 

As each UKIDSS frame has a different seeing value, it 
is necessary to break each fits file into its four constituent 
images. This is not necessary for SDSS images (which are 
stored in separate files). However, it is necessary to retrieve 
the SDSS image seeing from the image's tsField file. The 
SDSS image seeing is stored in the psf_width column of the 
tsField file. Where an image has a seeing better than our 
specified value, we use the fgauss utility to convolve our 
image down to our specified value. Where an image has a 
seeing worse than our specified value, we copy it without 
modification using the imcopy utility. Both utilities produce 
a set of UKIDSS files containing two HDUs: the original 
instrument header HDU and a single image HDU with seeing 
greater than or equal to our specified seeing. The output 
SDSS files contain just a single image HDU. This process 
takes approximated 2 seconds per frame. 



tion system. As UKIDSS images are stored in the ZPN 
projection system, SWARP internally converts the frames to 
the TAN projection system. There is also an astrometric 
distortion present in the UKIDSS images that SWARP 
corrects using the pv2-3 and pv2-l fits header parameter^]. 
SWARP is set to subtract the background from the image, us- 
ing a background mesh of 256 x 256 pixels (102 x 102 arcsec) 
and a back filter size of 3 x 3 to calculate the back- 
ground map. The backgroun d calculation follows the same 
algorithm as SExtractor (jBertin fc Arnoutsl 1 19961 '). To 
summarise: it is a bicubic-spline interpolation between the 
meshes, with a median filter applied to remove bright stars 
and artifacts. 

Every mosaic contains pixels that are covered by multiple 
input frames. SWARP is set to use the median pixel value 
when a number of images overlap. The effects of outlying 
pixel values, due to cosmic rays, bad pixels or CCD edges, 
should therefore be reduced. SWARP generates a weight map 
(Figure [4} that contains the flux variance in every pixel, 
calibrated using the background map described above. 
As the flux variance is affected by overlapping coverage, 
it is possible to see the survey footprint in the weight 
map. The weight map can be used within SExtractor to 
compensate for variations in noise. We do not use it when 
calculating our photometry for two reasons. Firstly, there is 
overlap between SDSS fpC frames. This overlap is not from 
observations, but from the method used to cut the long 
SDSS stripes into sections. SWARP would not account for 
this, and the weighting of the overlap regions on the optical 
mosaics would be calculated incorrectly. Secondly, using 
the weight maps would alter the effect of mosaic surface 
brightness limit variations upon our output catalogues. As 
we intend to model surface brightness effects later, we elect 
to use an unweighted photometric catalogue. 
A small number of objects will be split between input 
frames. SWARP can reconstruct them, with only small defects 
due to CCD edges. One such example is shown in Figure [5] 
We create both seeing-corrected and uncorrected mosaics 
for each passband and region combination. Each file is 20Gb 
in size. In total, the mosaics require just over 1 Terabyte of 
storage space. Each mosaic takes approximately 4 hours to 
create. 



3.4 Creation of master region images 

The SWARP utility is a multi-thread capable co-addition and 
image-warping too l that is part of the Terapix pipeline 
l|Bertin et alj|2002h . We use SWARP to generate complete im- 
ages of the GAMA regions from the normalised LAS/SDSS 
fits files. It is vital that the pixel size and area of coverage is 
the same for each filter, as SExtractor's dual-image mode 
requires perfectly matched frames. We define a pixel scale 
of 0.4 x 0.4 arcsec, and generate 117000 x 45000 pixel files 
centred around 09h00m00.0s, +01d00'00.0" (GAMA 9), 
12h00m00.0s, +00d00'00.0" (GAMA 12), and 14h30m00.0s, 
-l-OOdOO'00.0" (GAMA 15). SWARP is set to resample input 
frames using the default LANCZOS3 algorithm, which the 
Terapix team found was the optimal option when working 
with images from the Megacam instrument l|Bertin et al.l 
l2002h . 

SWARP produces mosaics that use the TAN WCS projec- 



4 PHOTOMETRY 

A major problem with constructing multi-wavelength cata- 
logues is that the definition of what constitutes an object can 
change across the wavelength range (see Appendix [A] par- 
ticularly Figure lA"Tj) . This can be due to internal structure 
such as dust lanes or star forming regions becoming brighter 
or fainter in different passbands, causing the extraction soft- 
ware to deblend an object into a number of smaller parts in 
one filter but not in another. This can lead to large errors in 
the resulting colours. We cannot be certain that the SDSS 
object extraction process would produce the same results 



1 An analysis of the astrometric distortion can be found in CASU 
doc ument VDF-TRE-IOA- 00009-0002 , currently available from 
: / / www . ast .cam .ac.uk/vdfs/docs/reports/ a strom / index.html | 



http 
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Figure 4. The r band weightmap of the 45000x45000 pixel subset 
region (5x5 deg; denned in Section [5}. Joins and overlap between 
frames are apparent (light grey). The mosaic does not have imag- 
ing of the top right corner or the bottom section (dark grey). 
These areas lie outside the region of interest as the mosaics are 
slightly larger in Declination than the GAMA regions themselves. 



Original UKIDSS image 



Master SWARP image 



Figure 5. A comparison between the original normalised image 
and the K band mosaic image of a galaxy on the bottom edge 
of an input UKIDSS frame. The bottom section of the galaxy is 
not part of this image and it has been stitched together on the 
mosaic using SWARP. 



as the extraction process we use to create our UKIDSS ob- 
ject catalogues. Seeing, deblending and aperture sizes will 
differ, compromising colours. To create a consistent multi- 
wavelength sample, the photometry needs to be recalculated 
consistently across all 9 filters. At the same time we can 
move from the circular apertures of SDSS and UKIDSS to 
full elliptical apertures, as well as investigate a variety of 
photometric methods. To gen erate our source catalog ues we 
use the SExtractor software (|Bertin &: Arnoutsll 19961 ). This 
is an object extraction utility, and its use is described in 
subsection 14.51 

In this paper we implement four different methods to de- 
fine our object positions and apertures. We produce three 



SExtractor catalogues (Ber tin fc Arnoutsl 19961) and on e 
Sersic catalogue (based upon GALFIT 3. IPeng et al.l I2007T ) . 
in addition to the original SDSS dataset. The generation of 
the three new SExtractor catalogues is detailed in section 
14.51 Each of the new SExtractor catalogues contain mag- 
nitudes calculated using two different elliptical, extended 
source apertures: the Kron and the Petrosian magnitude 
systems. They are briefly described in sections 14.11 14.21 and 
14.31 We also use a specially designed pipeline (SIGMA GAMA, 
iKelvin et al.ll201ol . based upon GALFIT 3) to calculate a total 
magnitude for each galaxy via its best fitting Sersic profile. 
This aperture system, and the process used to generate it, 
is described in Sections 14.41 and 14.61 

It is not obvious which photometric method will produce 
the optimal solution. Whilst the Sersic photometric method 
solves the missing light problem, it requires higher quality 
data to calculate the set of parameters that best model the 
galactic light profile. The Kron and Petrosian magnitude 
systems will work with lower quality data, but may under- 
estimate the flux produced by a galaxy. In this section we 
describe the photometric systems that we have used. Later, 
in sections [6] and [8] we will examine the different results 
produced by the choice of the photometric system. 



4.1 Self-defined Kron and Petrosian apertures 

We construct an independent catalogue for each filter, con- 
taining elliptical Kron and Petrosian apertures. These inde- 
pendent catalogues are then match ed across all 9 wavebands 
using STILTS (see section W7\ and iTavlorl l2006t ) . The aper- 
tures will therefore vary in size, potentially giving inconsis- 
tent colours, and as deblending decisions will also change, 
inconsistent matching between catalogues may occur. How- 
ever, as the apertures are defined from the image they are 
used on, there can be no problem with magnitudes being 
calculated for objects that do not exist, or with missing ob- 
jects that are not visible in the r band. 
The self-defined catalogues are generated from the basic mo- 
saics, where no attempt to define a common seeing for the 
mosaic has been made. This method should generate the 
optimal list of sources in each band; however, as the pre- 
cise definition of the source will vary with wavelength, the 
colours generated using this method will be inaccurate and 
subject to aperture bias. As the mosaic has variations in 
seeing, the PSF will also vary across the image. 



4.2 r band-defined Kron and Petrosian apertures 

We use SExtractor to define a sample of sources in the r 
band image, and then use the r band position and aperture 
information to calculate their luminosity within each filter 
(using the SExtractor dual image mode). As the aperture 
definitions do not vary between wavebands this method gives 
internally consistent colours, and as the list itself does not 
change source matching between filters is unnecessary. How- 
ever, where an object has changed in size (see Appendix |A"| . 
does not exist (e.g. an artifact in the r band sample) or when 
the r band aperture definition incorrectly includes multiple 
objects the output colours may be compromised. Any ob- 
ject that is too faint to be visible within the r band mosaic 
will also not be detected using this method. However, such 
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objects will be fainter than the GAMA sample's selection 
criteria, and would not be included within our sample. The 
r band-defined catalogues are generated from our seeing- 
degraded mosaics. They provide us with an optically-defined 
source sample. 

This method is analogous to the SDSS source catalogues, 
which define their apertures using the r passband data (un- 
less the object is not detected in r, in which case a dif- 
ferent filter is chosen). However, the GAMA photometric 
pipeline has a broader wavelength range as it now includes 
NIR measurements from the same aperture definition. Fur- 
thermore, the SDSS Petrosian magnitudes have not been 
seeing-standardised. While all data is taken at the same 
time, the diffraction limit is wavelength dependent and dif- 
ferent fractions of light will be missed despite the use of a 
fixed aperture. SDSS model magnitudes do account for the 
effects of the PSF. 

4.3 K band-defined Kron and Petrosian apertures 

This method works in the same way as the previous method, 
but uses the K band image as the detection frame rather 
than the r band image. We are limited in the total area, 
as the K band coverage is currently incomplete. However, 
for samples that require complete colour coverage in all 9 
filters, this is not a problem. As with r band-defined cata- 
logues, the K band-defined catalogues are generated from 
the seeing-corrected mosaics. They provide us with a NIR- 
defined source sample. The 7f-band defined Kron magni- 
tudes were used in the GAMA input catalogue (|Baldrv et al.l 
|2010T ) to calculate the star-galaxy separation J — K colour 
and the K band target selection. 

4.4 Sersic magnitudes 

We use the SIGMA modelling wrapper (see section 14.61 and 
iKelvin et ai]l20ld for more details) which in turn u ses the 
galaxy fitting software GALFIT 3.0 (|Peng et all 120071 1 to fit 
a single-Sersic component to each object independently in 
9 filters (ugrizY J H K), and recover their Sersic magni- 
tudes, indices, effective radii, position angles and elliptici- 
ties. Source positions are ini tially taken from the GAMA in- 
put catalogue, as defined in iBaldrv et all (|2010h . All Sersic 
magnitudes are self-defined; as each band is modelled inde- 
pendently of the others, the aperture definition will vary and 
colour may therefore be compromised. 

Single-Sersic fitting is comparable to the SDSS model mag- 
nitudes. SIGMA therefore should recover total fluxes for ob- 
jects that have a Sersic index in the range 0.3 < n < 20, 
where model magnitudes force a fit to either an exponen- 
tial (n=I) or deVaucouleurs (n=4) profile. The systematic 
magnitude errors that arise when model magnitudes are fit 
to galaxies that do not follow an exponential or deVau- 
couleurs profile ijGrahaml [ioOll : H rown et al.ll2003h do not 
occur in SIGMA. The SDSS team developed a composite mag- 
nitude system, CMODEL, that calculates a magnitude from 
the combination of the n=l and n=4 syste ms, in order to cir- 
cumvent this issue ()Abazaiian et al1l2004l ). We compare our 
Sersic magnitudes to their results later. Sersic magnitudes do 
not suffer from the missing-flux issue that affects Kron and 
Petrosian apertures. Petrosian magnit udes may underest i- 
mate a galaxy's luminosity by 0.2 mag l|Strauss et al.i r2002). 



while under certain conditions a Kron a perture may on ly re- 
cover half of a galaxy's total luminosity (I Andreonl2 002) . The 
Sersic catalogues are generated from the seeing-uncorrected 
mosaics, as the seeing parameters are modelled within SIGMA 
using the PSFEx software utility (E. Bertin, priv. comm). 



4.5 Object Extraction of Kron and Petrosian 
apertures 

The SExtractor utility (Ber tin fc Arnoutsi|l996l ) is a pro- 
gram that generates catalogues of source positions and aper- 
ture fluxes from an input image. It has the capacity to 
define the sources and apertures in one frame and calcu- 
late the corresponding fluxes in a second frame. This dual 
image mode is computationally more intensive than the 
standard SExtractor single image mode (in single image 
mode, SExtractor can extract a catalogue from a mosaic 
within a few hours; dual image mode takes a few days 
per mosaic). Using the u, g, r, i, z, Y, J, H and K im- 
ages created by the SWARP utility, we define our catalogue 
of sources independently (for the self-defined catalogues), 
using the r band mosaics (for the r band-defined cata- 
logue) or the K band mosaics (for the K band-defined cat- 
alogue) and calculate their flux in all nine bands. The nor- 
malisation and SWARP processes removed the image back- 
ground and standardised the zeropoint; we therefore use 
a constant MAG.ZEROPOINT=30, and BACK.VALUE=0. 
SExtractor generates both elliptical Petrosian (2.0 -Rpetro) 
and Kron- like apertures (2.5 Rktou, called AUTO magni- 
tudes). SExtractor Petrosian magnitudes are computed us- 
ing — - — = 0.2, the same parameter as SDSS. As the 
mosaics have been transformed onto the AB magnitude sys- 
tem, all magnitudes generated by the GAMA photometric 
pipeline are AB magnitudes. 

The seeing convolution routine smooths out the background 
and correlates the read noise of the images (this is apparent 
in Figure. lATj) . As SExtractor detects objects of > xa above 
the background (where a; is a definable parameter, set to 1 in 
the default file and for our seeing unconvolved catalogues), 
this assists the detection process, allowing SExtractor to 
find objects to a much greater depth, thus increasing the 
number of sources extracted using the standard setup. How- 
ever, these new objects are generally much fainter than the 
photometric limits of the GAMA spectroscopic campaign, 
many are false detections, and the time required to generate 
the source catalogues (particularly using SExtractor in dual 
image mode) is prohibitively large. Using a 10000x10000 
pixel subset of the GAMA9 r band mosaic, we have at- 
tempted to calculate the DETECT. THRESH parameter 
that would output a catalogue of approximately the same 
depth and size as the unconvolved catalogue within our spec- 
troscopic limits (see section [2~Tj) . The distribution of objects 
with different DETECT_THRESH sigma parameters, com- 
pared to the unconvolved catalogue, is shown in Figure[6] We 
use a DETECT. MIN AREA of 9 pixels. As the unconvolved 
catalogue is slightly deeper than the 2a, but not as deep as 
the 1.7a convolved catalogue, we use a DE TEC T_ THRESH 
parameter of 1 .7a to generate our convolved catalogues. The 
1.7a and la catalogues have consistent number counts to 
?"auto = 21 mag; half a magnitude beyond the tsdss band 
magnitude limit of the GAMA input catalogue. 
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Figure 6. The effects of changing the SExtractor DE- 
TECT_THRESH parameter on a subset of an r band mosaic. 
The dotted black line is the deepest r band sample limit of the 
GAMA survey (for those objects that are K or z selected). 



4.6 Object extraction for Sersic magnitudes 

Sersic magnitudes are obtained as an output from the 
galaxy modelling program SIGMA (Structural Investigation 
of Galaxies via Model Analysis) wr itten in the R program- 
ming language (|Kelvin et al.l l20l"oh . In brief, SIGMA takes 
the Right Ascension and Declination of an object that has 
passed our star-galaxy separation criteria and calculates its 
pixel position within the appropriate mosaic. A square re- 
gion, centred on the object, is cut out from the mosaic con- 
taining a minimum of 20 guide stars with which to generate 
a PSF. SExtractor then provides a FITS-standard input 
catalogue to PSF Extractor (E. Bertin, priv. comm.) which 
generates an empirical PSF for each image. Ellipticities and 
position angles are obtained from the STSDAS Ellipse 
routine within IRAF, and provides an input to Galf it. The 
larger cutout is again cut down to a region which contains 
90% of the target object's flux plus a buffer of 10 pixels, 
and will only deviate from this size if a bright nearby object 
causes the fitting region to be expanded in order to model 
any satellites the target may have. 

GALFIT 3 is then used to fit a single-Sersic component to 
each target, and several runs may be attempted if, for exam- 
ple, the previous run crashed, the code reached its maximum 
number of iterations, the centre has migrated to fit a sep- 
arate object, the effective radius is too high or low or the 
Sersic index is too large. SIGMA employs a complex event 
handler in order to run the code as many times as neces- 
sary to fix these problems, however not all problems can be 
fixed, and so residual quality flags remain to reflect the qual- 
ity of the final fit. The SIGMA package takes approximately 
10 seconds p er object. For f ull de tails of the SIGMA modelling 
program, see iKelvin et all (|2010h . 



4.7 Catalogue matching 

The definition of the GAMA spectroscopic target selection 
(herein refer r ed to as the tiling catalogue) is detailed in 
iBaldrv et all (|2010h . and is based on original SDSS DR6 
data. We therefore need to relate our revised photometry 
back to this catalogue in order to connect it to our AAOmega 
spectra. The tiling catalogue utilises a mask around bright 
stars that should remove most objects with bad photometry 
and erroneously bright magnitudes, as well as implement- 
ing a revised star-galaxy separation quantified against our 
spectroscopic results. It has been extensively tested, with 
sources that are likely to be artifacts, bad deblends or sec- 
tions of larger galaxies viewed a number of times by different 
people. By matching our catalogues to the tiling catalogue, 
we can access the results of this rigorous filtering process, 
and generate a full, self-consistent set of colours for all of the 
objects that are within the GAMA sample (and are within 
regions that have been observed in all nine passbands). As 
the tiling catalogue is also used when redshift targeting, we 
will be able to calculate the completeness in all the pass- 
bands of the GAMA survey. The GAMA tiling catalogue is 
a subset of the GAMA master catalogue (herein referred to 
as the master catalogue). The master catalogue is created 
using the SDSS DR6 catalogue stored on the CAlfl Un- 
like the master catalogue, the tiling catalogue undertakes 
star-galaxy separation, and applies surface brightness and 
magnitu de selection. 

STILTS (|Tavlorll2006l ) is a catalogue combination tool, with 
a number of different modes. We use it to join our region 
catalogues together to create r-defined, A"-defined and self- 
defined aperture photometry catalogues that cover the entire 
GAMA area. We also use it to match these catalogues to the 
GAMA tiling catalogue. 



4.8 Source catalogues 

The catalogues that have been generated are listed in 
Table [2] The syntax of the Key column is as follows. 
X[u] means a u band magnitude from an X band-defined 
aperture, {u} means a self-defined u band magnitude 
and + denotes a STILTS tskymatch2 5 arcsec, unique 
nearest-object match between two catalogues (see Section 
I4.7|l . Where two datasets are combined together without 
the + notation (i.e., the final two lines), this denotes 
a STILTS tmatch2, matcher = "exact" match using SDSS 
OBJID as the primary key. Note that in a set of self-defined 
samples ({ugrizY J H K}) , each sample must be matched 
separately (as each contains a different set of sources), 
and then combined. This is not the case in the aperture 
defined samples (where each sample contains the same set 
of sources). Subscripts denote the photometric method used 



2 We use the query SELECT * FROM dr6.PhotoObj as p 
WHERE ( p.modelmag_r - p. extinctions < 20.5 or p.petromag_r 
- p. extinctions < 19.8 ) and ( (p.ra > 129.0 and p.ra < 141.0 and 
p. dec > -1.0 and p.dec < 3.0) or (p.ra > 174.0 and p.ra < 186.0 
and p.dec > -2.0 and p.dec < 2.0) or (p.ra > 211.5 and p.ra < 
223.5 and p.dec > -2.0 and p.dec < 2.0) ) and ((p. mode = 1) or 
(p. mode = 2 and p.ra < 139.939 and p.dec < -0.5 and (p. status 
& dbo.fphotostatus('OK.SCANLINE')) > 0)) 
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Catalogue Name Key Abbreviation 

r-defincd catalogue r[ugriz]sDSS + r[ugrizY J H K\qam A-.Petro,Kron catrdef 

self-defined catalogue r[ugriz]sDSS + {ugrizY J H K}Q A MA:Petro,Kron catsd 

/^-defined catalogue r[ugriz\sDSS + K[ugrizY J H K]q A m A-.Pe.tro Kron catKdef 

Sersic catalogue r[ugriz]sDSs{ u 9 r "i- z YJHK}c A MA:Seraic catsers 
GAMA master catalogue (r[ugriz} SDSS {r} GAMA . Sl , rsic ) + {ugrizY J H K} GAMA . Petro 

,Kron catmast 

Table 2. The names of the generated catalogues, the prescription used to create them and their abbreviated filename. The syntax in 
the Key column is summarised in section [4,81 



for each catalogue. 



5 TESTING THE GAMA CATALOGUES 

In order to test the detection and deblending outcomes 
within the GAMA catalogues, a subsection of 25 sq deg has 
been chosen from near the centre of the GAMA 9 region (the 
pixels used are 20000-65000 in the x direction of the mosaic, 
and 0-45000 in the y direction). This region was chosen as 
it contains some of the issues facing the entire GAMA sub- 
set, such as area incompleteness. UKIDSS observations miss 
a large fraction of the subset area - approximately 3.02 sq 
deg of the region has incomplete NIR coverage. The subset 
region was also chosen because it partially contains area cov- 
ered by the Hersc hel ATLAS science verification region (see 
lEales et al J 120091 ) . Within this region, we ran SExtractor, 
and compared our results with the source lists produced by 
the SDSS and UKIDSS extraction software. Unless other- 
wise specified, all magnitudes within this section were cal- 
culated using r-defined apertures. 

5.1 Numerical breakdown 

After generating source catalogues containing self-consistent 
colours for all objects in the subset region (using the process 
described in subsections 14.51 and 14. Tf> . we are left with an r 
band aperture-defined subset region catalogue containing 
1810134 sources and a K band aperture-defined subset 
region catalogue containing 2298224 sources (these are 
hereafter referred to as the r band and K band catalogues). 
These catalogues contain many sources we are not interested 
in, such as sources with incomplete colour information, 
sources that are artifacts within the mosaics (satellite trails, 
diffraction spikes, etc), sources that are stars and sources 
that are fainter than our survey limits. 

The unfiltered r and K band catalogues were matched to 
the master catalogue with a centroid tolerance of 5 arcsec, 
using the STILTS tskymatch2 mode (see Section I4.7[l . 
Table [3] contains a breakdown of the fraction of matched 
sources that have credible X A uto and Xpetro for all 
nine passbands (sources with incorrect AUTO or PETRO 
magnitudes have the value 99 as a placeholder; we impose 
a cut at X = 50 to remove such objects). Generally, the 
low quality of the u band SDSS images causes problems 
with calculating extended source magnitudes, and this 
shows itself in the relatively high fraction of incomplete 
sources. This problem does not affect the other SDSS filters 
to anywhere near the same extent. SDSS observations 
do not cover the complete subset area, but they have 
nearly complete coverage in both the r-defined (which is 



Band % Cover Sources (r) % (r) Sources (K) % (K) 



Total 




129488 




123740 




u 


100 


111403 


86 


105801 


86 


g 


100 


129169 


100 


123317 


100 


r 


100 


129481 


100 


123610 


100 


i 


100 


129358 


100 


123533 


100 


z 


100 


128287 


99 


122479 


99 


Y 


88 


108167 


84 


109672 


89 


J 


89 


109364 


84 


109816 


89 


11 


96 


121212 


94 


121846 


98 


K 


94 


118224 


91 


122635 


99 



Table 3. Number of sources within the subset region with good 
SExtractor X Au t and Xp e t ro , where X is ugrizY J H K , from 
the r or K band-defined aperture catalogues matched to the 
GAMA master catalogue. The total number of sources within 
the GAMA master catalogue for this region of sky is 138233. % 
Cover is defined relative to r band cover; where SDSS coverage 
does not exist there are no GAMA master catalogue sources. 

dependent on SDSS imaging) and -ft'-defined (reliant on 
UKIDSS coverage) catalogues. The UKIDSS observations 
cover a smaller section of the subset region, with the Y 
and J observations (taken separately to the H and K) 
covering the least area of sky. This is apparent in the 
r band catalogue, where at least 16% of sources lack 
PETRO or AUTO magnitudes in one or more passband. 
By its definition the K band catalogue requires K band 
observations to be present; as such there is a high level of 
completeness in the grizH and K passbands. However, the 
number of matched SDSS sources in the K band catalogue 
itself is 4.2% lower than in the r band catalogue. 

There are 138233 master catalogue SDSS sources within 
the subset region. 119330 SDSS objects have matches 
(within a 5 arcsec tolerance) in both the r band and K 
band master-cat matched catalogues (this number is found 
by matching SDSS OBJID between the catalogues). Those 
SDSS objects that do not have matches in both master-cat 
matched catalogues are shown in Figure [7] We detail the 
reasons for the missing objects in section \57Z\ 



5.2 SDSS sources missing in the master 

catalogue- aperture matched catalogues 

There are 18903 SDSS sources that are not found when the 
master catalogue is matched to either the r or K-de&ned 
subset region catalogues; 13.7% of the total number of 
master catalogue sources within the subset region. Figure [7] 
shows their distribution on the sky. 8745 sources are not 
found within the r-defined catalogue (6.3% of the master 
catalogue sample) and 14493 are not found within the 
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Figure 7. The Right Ascension and Declination of SDSS ob- 
jects that are not in either the r or K band master-cat matched 
catalogue. The darker areas notes a high density of unmatched 
objects. 



Reason for non-detection 


Number 
of objects 


% of GAMA master cata- 
logue subset region sample 


Possible dcblcnding mis- 


601 


0.4 


match 






Saturation spike / satel- 


404 


0.3 


lite or asteroid trail 






SExtractor detection 


3831 


2.8 


failure 






Either a low surface 


2391 


1.7 


brightness source or no 






source 






Part of a large deblcnded 


1463 


1.1 


source 






Low image quality mak- 


55 


0.04 


ing detection difficult 







Table 4. A breakdown of the reasons for faulty detections in the 
8745 SDSS objects that are not matched to the r band subset 
region catalogue. The images of the SDSS objects were generated 
from the standard r band GAMA mosaics, and all 8745 objects 
were viewed by one observer (DTH). The criteria selection is as 
follows. The first category is chosen in those cases where an object 
has a nearby neighbour or may have been deblended into multiple 
sources by the SDSS algorithm. The second category is chosen 
where the position of the object is covered by a spike/trail. The 
third category is where a source is visible by eye. The fourth 
category is where a source is not visible above the noise. The fifth 
category is chosen when a source is obviously part of a larger 
structure. The sixth category is chosen when the SDSS data is 
too low quality for visual classification to be undertaken. 



if-defined catalogue (10.5% of the sample), with 4335 of 
the sources unmatched to either the r or the /^-defined 
sample (3.1% of the master catalogue sample). As the 
SDSS sample is defined by optical data, it is unsurprising 
that a far larger number of sources are not found within 
the if-defined catalogue. Of the 18903 unmatched master 
catalogue sources, only 2367 have passed star-galaxy 
separation and are brighter than the GAMA spectroscopic 
survey magnitude limits (r < 19.8 or zK selected). 

Using r band imaging, we have visually inspected all 
8745 SDSS sources where our r-defined subset region cata- 
logue cannot find a match within 5 arcsec. Table [4] contains 
a summary of the reasons we do not find a r band match. 
Using the SExtractor detection failure rate from the sub- 
set region as a guide to the detection failure rate for the 
entire GAMA region, SExtractor will miss approximately 
2.8% of the objects recovered by SDSS. A second problem 
was flagged through the inspection process; a further 1.7% of 
the master catalogue sources within the subset region were 
not visible. Either these objects are low surface brightness 
extended objects, possibly detected in a different band, or 
the SDSS object extraction algorithm has made a mistake. 
A further 1.8% of sources within the GAMA master cat- 
alogue will be missed by SExtractor due to differences in 
deblending decisions (either failing to split two sources or 
splitting one large object into a number of smaller parts), 
low SDSS image quality making SExtractor fail to detect 
any objects, or an artifact in the image being accounted for 
by SExtractor (such as a saturation spike from a large star 
being detected as a separate object in SDSS). 



5.3 Sources in our r band catalogue that are not 
in the GAMA master catalogue 

To be certain that the SDSS extraction software is giving 
us a complete sample, we check whether our r band sub- 
set region catalogue contains sources that should be within 
the GAMA master catalogue but are not. There are 61351 
sources within the r-defined subset region catalogue that 
have a complete set of credible AUTO and PETRO magni- 
tudes, and are brighter than the GAMA spectroscopic sur- 
vey limits. 619 of these sources do not have SDSS coun- 
terparts. We have visually inspected these sources; a break- 
down is shown in TableJS] Similar issues cause missing detec- 
tions using the SDSS or SExtractor algorithms. However, 
some of the unseen sources that SExtractor detected may 
be due to the image convolution process (Section 13. 3p gath- 
ering up the flux from a region with high background noise 
and rearranging it so that it overcomes the detection thresh- 
old. Figure [5] shows the distribution of rAuto ^ 20.5 mag 
sources detected when SExtractor is run upon an original 
SDSS image file (covering ~ 0.04 deg 2 ), and the sources from 
the same file after it has undergone the image convolution 
process. 233 sources are found in the original SDSS frame, 
and 3 additional sources are included within the convolved 
frame sample. An examination of two sources that are in 
the convolved frame dataset and not in the original sam- 
ple shows the effect: these sources have rAuto luminosities 
of 20.64 mag and 20.77 mag pre-convolution, but rAuto lu- 
minosities of 20.40 mag and 20.48 mag post-convolution. 
Taking the SDSS non-detection rate within the subset re- 
gion to be the same as the non-detection rate over the entire 
GAMA region, we expect that the SDSS algorithm will have 
failed to detect 0.1% of sources brighter than the GAMA 
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Type of source Number of objects 
Source 171 
No visible source 274 
Section of bright star 163 
Possible deblend mismatch 10 
Low image quality making detection difficult 1 



6 PROPERTIES OF THE CATALOGUES 
6.1 Constructing a clean sample 

In order to investigate the photometric offsets between 
different photometric systems, we require a sample of 
galaxies with a complete set of credible photometry that are 
unaffected by deblending decisions. This has been created 
via the following prescription. We match the r-defined 
aperture catalogue to the GAMA master catalogue with a 
tolerance of 5 arcsec. We remove any GAMA objects that 
have not been matched, or have been matched to multiple 
objects within that tolerance (when run in All match mode 
STILTS produces a GroupSize column, where a NULL value 
signifies no group). We then match to the 9 self-defined 
object catalogues, in each case removing all unmatched 
and multiply matched GAMA objects. As our convolution 
routine will cause problems with those galaxies that contain 
saturated pixels, we also remove those galaxies that are 
flagged as saturated by SDSS. This sample is then linked 
to the Sersic pipeline catalogue (using the SDSS objid as 
the primary key). We remove all those Sersic magnitudes 
where the pipeline has flagged that the model is badly 
fit or where the photometry has been compromised and 
match to the K band aperture-defined catalogue, again 
with unmatched and multiple matched sources removed. 
This gives us a final population of 18065 galaxies that 
have clean r-defined, ^-defined, self-defined and Sersic 
magnitudes, are not saturated and cannot have been 
mismatched. Having constructed a clean, unambiguous 
sample of common objects, any photometric offset can only 
be due to differences between the photometric systems 
used. As we remove objects that are badly fit by the Sersic 
pipeline, it should be noted that the resulting sample will, 
by its definition, only contain sources that have a light 
profile that can be fitted using the Sersic function. 



Table 5. A breakdown of the 619 r-defined subset region cat- 
alogue objects brighter than the GAMA sample limits that are 
not matched to the GAMA master catalogue. The images of the 
subset region catalogue objects were generated from the standard 
r band GAMA mosaics, and all 619 objects were viewed by one 
observer (DTH). 




Figure 8. A comparison between the objects detected when 
SExtractor is run over an original SDSS image, and when it is 
run over the convolved, mosaic imaging. Yellow circles are sources 
with r 20.5 mag detected from the GAMA mosaic, red crosses 
are sources that are detected from the original SDSS data. 

spectroscopic survey limits; approximately 1000 sources will 
not have been included within the master catalogue. 



5.4 Sources in our K band catalogue that are not 
in the UKIDSS DR5PLUS database 

We have also tested the UKIDSS DR5 catalogue. We 
have generated a catalogue from the WSA that selects all 
UKIDSS objects within the GAMA subset regior0, and we 
have matched this catalogue to the K band-defined subset 
region catalogue. From the 69537 K band-defined subset 
region catalogue sources, there are 4548 sources that have 
not been matched to an UKIDSS object within a tolerance 
of 5 arcsec. We have visually inspected K band images of 
those objects that are brighter than the GAMA spectro- 
scopic survey K band limit (Kauto ^ 17.6 mag). We find 
that 29 of the 117 unmatched objects are real sources that 
are missed by the UKIDSS extraction software; a negligi- 
ble fraction of the entire dataset. A large (but unquantified) 
fraction of the other 88 sources are suffering from the convo- 
lution flux-redistribution problem discussed in Section 15.31 
The background fluctuations in K band data are greater 
than in the r band, making this a much greater problem. 

3 We use the query "SELECT las.ra, las. dec, las.kPetroMag 
FROM lasSource as las WHERE las.ra < 139.28 AND las.ra > 
134.275 AND las.dec > -1 AND las.dec < 3" 



6.2 Photometric offset between systems 



Figures [9l 1 101 [TT1 1 121 131 and 1 141 show the dispersion between 
different photometric systems produced by this sample. In 
Figure [5] we compare Kron and Petrosian magnitudes; in all 
other figures we compare the photometric system to SDSS 
petromag. In all photometry systems, the gri relationships 
are tightest, with the u and z relationships subject to a 
greater scatter, breaking down almost entirely for Figures 
[T2l and \T3\ The correlation between the SDSS petromag 
and the r-defined Petrosian magnitude (Figure lll[) looks 
much tighter than that between the SDSS petromag and 
the self-defined Petrosian magnitude (Figure [T^J. The stan- 
dard deviation of the samples are similar, with marginally 
more scatter in the self-defined sample (0.129 mag against 
0.148 mag). The median offset between SDSS Petrosian and 
the r-defined Petrosian magnitude, however, is 0.01 mag 
greater. 

Figure 1131 illustrating the relationship between the Sersic 
magnitude and the SDSS petromag, produces median 
AmsDSS — msersic values of 0.12, 0.06, 0.06, 0.07 and 
0.09 mag in ugriz. Thes e values can be comp ared to those 
presented in Figure 13 of lBlanton et all i|2003T ) (-0.14, 0.00, 
0.06, 0.09 and 0.14 mag at z — 0.1, using the 01 ugriz filters), 
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Figure 9. GAMA r-defined aperture Petrosian minus Auto magnitudes for a clean sample of galaxies in ugriz. Contours are for 4 to 
512 galaxies per bin, rising geometrically in powers of 2. Bins are 0.05 magx 0.05 mag in size. 
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Figure 10. SDSS petromag minus GAMA r-defined aperture Auto magnitudes for a clean sample of galaxies in ugriz. Contours are 
for 4 to 512 galaxies per bin, rising geometrically in powers of 2. Bins are 0.05 magx0.05 mag in size. 
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Figure 11. SDSS petromag minus GAMA r-defined aperture Petrosian magnitudes for a clean sample of galaxies in ugriz. Contours 
are for 4 to 512 galaxies per bin, rising geometrically in powers of 2. Bins are 0.05 magx0.05 mag in size. 
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Figure 12. SDSS PETROMAG minus GAMA self-defined aperture Petrosian magnitudes for a clean sample of galaxies in ugriz. Contours 
are for 4 to 512 galaxies per bin, rising geometrically in powers of 2. Bins are 0.05 magx0.05 mag in size. 
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Figure 13. SDSS petromag minus GAMA Sersic magnitudes for a clean sample of galaxies in ugriz. Contours are for 4 to 512 galaxies 
per bin, rising geometrically in powers of 2. Bins are 0.05 magx0.05 mag in size. 
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Figure 14. SDSS CMODEL minus GAMA Sersic magnitudes for a clean sample of galaxies in ugriz. Contours are for 4 to 512 galaxies 
per bin, rising geometrically in powers of 2. Bins are 0.05 magx0.05 mag in size. 
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given the variance in the relationship (the standard devia- 
tion in our samples are 0.77, 0.28, 0.21, 0.22 and 0.40 mag, 
respectively). A significant fraction (~ 28%) of the sam- 
ple has rsDSS — rsersicl.0.5 mag, and therefore lies beyond 
the boundaries of this image. These offsets are significant, 
and will be discussed further in Section [7] We can say that 
the r-defined aperture photometry is the closest match to 
SDSS petromag photometry. Figure [14] shows the relation- 
ship between the GAMA Sersic magnitude, and the optimal 
model magnitude provided by SDSS (cmodel). The model 
magnitudes match closely, with negligible systematic offset 
between the photometric systems in gri. 

6.3 Colour distributions 

In order to identify the optimal photometric system, we as- 
sume that intrinsic colour distribution of a population of 
galaxies can be approximated by a double-Gaussian distri- 
bution (the superposition of a pair of Gaussian distributions 
with different mean and standard deviation parameters). 
This distribution can model the bimodality of the galaxy 
population. The presence of noise will broaden the distri- 
bution; hence the narrowest colour distribution reveals the 
optimal photometric system for calculating the colours of 
galaxies, and therefore deriving accurate SEDs. Figure [15] 
shows the (u — r) and (r — K) colour distributions for each 
photometric system, for objects within our subset region. 
In order to calculate the dispersion in the colour distribu- 
tion, we generate colour-distribution histograms (with bins 
of 0.1 mag), and find the double-Gaussian distribution pa- 
rameters that best fit each photometric system. The best- 
fitting standard deviation parameters for each sample are 
shown at the bottom of each plot, and are denoted ax,i and 
o~x,2 (where X is the photometric system fitted). The sam- 
ple with the smallest set of a parameters should provide the 
optimal photometric system. 

The SDSS, GAMA r-defined aperture and GAMA K- 
defined distributions (the first, third and fourth diagrams on 
the top two rows) show a very similar pattern; a tight dis- 
tribution of objects with a small number of red outliers. As 
expected, when we use apertures that are defined separately 
in each filter (the second diagram on the top two rows), 
the colour distribution of the population is more scattered 
[crpetro,! = 0.7576 mag, <Tp etro ,2 = 0.7919 mag, a A uto,i = 
0.5886 mag, <7Auto,2 = 0.7086 mag) and does not show the 
bimodality visible in the matched aperture photometry (at 
the bright end of the distribution there are two distinct sub- 
populations; one sub-population above u — r = 2 mag, the 
other below). For the same reason, and probably because 
of the low quality of the observations, the (u — r) plot us- 
ing the Sersic magnitudes (the final diagram on the top row) 
has the broadest colour distribution (crsersic.l = 0.6242 mag, 
o~Sersic,2 ~ 1.098 mag), although it is well behaved in (r-K). 
To generate a series of (r — K) colours using the UKIDSS 
survey (leftmost plot on the bottom two rows), we have 
taken all galaxies within the UKIDSS catalogujfl and match 

4 We run a query at the WSA on UKIDSSDR5PLUS looking for 
all objects within our subset region with lasSource.pGalaxy > 
0.9 & lasSource.kPetroMag < 20 - equivalent to Kab < 
21.9 mag 



them (with a maximum tolerance of 5 arcsec) to a copy of 
the tiling catalogue that had previously been matched with 
the K band aperture-defined catalogue. The distribution of 
(r — K) colours taken from the SDSS and UKIDSS survey 
catalogues is the first diagram on the bottom two rows of 
the image. As the apertures used to define the UKIDSS and 
SDSS sources are not consistent, we find that the tightest 
(r — K) distribution comes from the GAMA Jf-defined aper- 
ture sample (fourth from the left on the bottom row, with 
(JAuto,i = 0.3137 mag, a A uto,2 = 0.4921 mag). The GAMA 
sample that relies on matching objects between self-defined 
object catalogues (the second diagram on the bottom two 
rows) has the broadest distribution (o\p etr0i i = 0.3359 mag, 
op e tro,2 = 0.6015 mag). The distribution of sources in the 
Sersic (r — K) colour plot is much tighter than in the 
(u — r), though still not as tight as the distribution in the 
fixed aperture photometric systems (asersic.i = 0.364 mag, 
<TSersic,2 = 0.6159 mag). Figure [151 confirms the utility of 
the GAMA method: by redoing the object extraction our- 
selves, we have generated self-consistent colour distribu- 
tions based on data taken by multiple instruments that 
has a far smaller scatter than a match between the sur- 
vey source catalogues ((tsds.s+ukidss,i = 0.3342 mag, 
&sdss+ukidss,2 = 0.5807 mag). 

We provide one more comparison between our colour dis- 
tribution and that provided by SDSS and UKIDSS survey 
data. Figure [16] displays the X — H distribution produced by 
the GAMA galaxies with complete ugrizY J H K photome- 
try and good quality redshifts within 0.033 < z < 0.6. The 
effective wavelengths of the filter set for each galaxy are 
shifted using the redshift of the galaxy. The colour distri- 
bution provided by the GAMA photometry produces fewer 
outliers than the SDSS /UKIDSS survey data sam ple, and is 
well constrained by the iBruzual fc Charlotl ||2003h models. 



7 FINAL GAMA PHOTOMETRY 

Sections [5] and [6] show that the optimal deblending out- 
come is produced by the original SDSS data, but the best 
colours come from our r-defined aperture photometry (Sec- 
tion 16. 3[) . We see that our r-defined aperture photometry 
agrees with the SDSS petromag photometry. However, we 
have also demonstrated that SDSS petromag misses flux 
when compared to our Sersic total magnitude. Here we com- 
bine these datasets to arrive at our final photometry. We 
combine the SDSS deblending outcome with our r-defined 
aperture colours and the Sersic total-magnitude to produce 
our best photometric solution. 

7.1 Sersic magnitudes 

To check the reliability of the Sersic photometry pipeline, 
we must examine its distribution against a photometric sys- 
tem we consider reliable. We examine the distribution of the 
Sersic photometry against our r-defined AUTO photometry. 
Figure IT7I shows the distribution of Sersic - GAMA r-defined 
AUTO magnitude against r-defined AUTO magnitude for 
all objects in the GAMA sample that have passed our star- 
galaxy separation criteria and have credible AUTO mag- 
nitudes. Whilst there is generally a tight distribution, the 
scatter in the u band, in particular, is a cause for concern. 
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Figure 15. A comparison between the u minus r and r minus K colours produced using SDSS MODELMAG, GAMA self-defined Pctrosian 
magnitudes, GAMA r / K-defined Petrosian magnitudes, Sersic and GAMA Total magnitudes for objects in the subset region. Contours 
shown are 2 to 512 galaxies per bin, rising geometrically in powers of 2. Bins are 0.1 mag in width in each axis. The a parameter comes 
from the best-fitting bivariate-Gaussian distribution, when it is fit to the colour-distribution histogram in each plot. 



GAMA: The photometric pipeline 21 





Figure 16. A comparison between the X — H colours produced using SDSS magnitudes and GAMA Total magnitudes. Data comes 
from all GAMA galaxies with good quality redshifts (0.033 < z < 0.6) and complete ugrizY J H K photometry. Effective wavelengths are 
calculated from the redshift of the galaxy and the filter effective wavelength, and the datas et is binned int o a 50x50 bin matrix. Two 
Bruzual-Charlot 03 SSP instantaneous-burst models are also plotted. Both models use the IChabrierl l|2003l l IMF, with mass cutoffs at 
O.f and fOOMQ. Stellar evolution is undertaken using the Padova f994 prescription. The dark grey line is a model evolved to If Gyr , 
with Z=0.05 and Y=0.352. The purple line is a model evolved to 0.25 Gyr, using Z=0.02 (Z©) and Y=0.28. 



iGraham et ail |2005h analytically calculate how the ratio of 
Sersic flux to Petrosian flux changes with the Sersic index 
of the object. The fraction of light missed by a Petrosian 
aperture is dependent upon the light profile of source. Fig- 
ure [19] shows the distribution of Sersic - GAMA r-defined 
Petrosian magnitude against Sersic index, redshift, absolute 
and apparent magnitude for all r-band objects in the GAMA 
sample that have passed our star-galaxy separation criteria, 
and have credi ble r, u and K r-defined PETRO magnitudes. 
IGraham et al.l report a 0.20 mag offset for an n = 4 pro- 
file, and a 0.50 mag offset for an n = 8 profile. The median 
rsersic — rp etr0 sian offset for objects with 3.9 < n < 4.1 in 
this sample is —0.115 mag, with rms scatter of 0.212 mag, 
and —0.408 mag, for objects with 7.9 < n < 8.1, with rms 
scatter of 0.292 mag. Both results agree with the reported 



values. We have plotted the magnitude offset with Sersic in- 
dex fu nction from Figure. 2 (their Panel a) of lGraham et al.l 
(2005) in the uppermost plot of Figure 1191 The function 
is an extremely good match to our photometry. Figure 1181 
shows the distribution of Sersic - SDSS CMODEL magnitude 
against Sersic index, redshift, absolute and apparent magni- 
tude for all r-band objects in the GAMA sample that have 
passed our star-galaxy separation criteria, and have credible 
r, u and K r-defined PETRO magnitudes. The distributions 
are very similar to those produced by the Sersic-Petrosian 
colours in Figure [191 An exception is the distribution with 
Sersic index, where the Sersic - CMODEL offset is distributed 
closer to mag, until n=4, at which point the Sersic magni- 
tude detects more flux. As the CMODEL magnitude is defined 
as a combination of n=l and n=4 profiles, it is unsurpris- 
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ing that it cannot model high n profile sources as well as 
the GAMA Sersic magnitude, which allows the n parameter 
greater freedom. 

The r band Sersic magnitude shows no anomalous be- 
haviour. Sersic profiling is reliable when undertaken using 
the higher quality SDSS imaging (particularly gri), but not 
when using the noisier u band data. It is clear that the 
u band Sersic magnitude is not robust enough to support 
detailed scientific investigations. In order to access a Sersic- 
style total magnitude in the u band, we are therefore forced 
to create one from existing, reliable data. We devise such an 
approach in Section [7721 



7.2 'Fixed aperture' Sersic magnitudes 



As mentioned in Section 14.41 the Sersic magnitude is taken 
from a different aperture in each band. We therefore cannot 
use Sersic magnitudes to generate accurate colours (compare 
the scatter in the Sersic colours and the AUTO colours in 
Figure [20]) . We also do not consider the u band Sersic mag- 
nitudes to be credible (see Section I7.1[l . However, we also 
believe that the r band Sersic luminosity function may be 
more desirable than the light-distribution defined aperture 
r band luminosity functions. The calculation of the total 
luminosity density using a non-Sersic aperture system may 
underestimate the parameter. We require a system that ac- 
counts for the additional light found by the Sersic magni- 
tude, but also provides a credible set of colours. 
We derive a further magnitude X to tai, using the equation 
Xtotai = (Xauto — r a uta) + r Sersic, where auto is the r- 
defined AUTO magnitude. In effect, this creates a measure 
that combines the total r band flux with optimal colours, us- 
ing SDSS deblending to give us the most accurate catalogue 
of sources (by matching to the GAMA master catalogue); 
the best of all possibilities. We accept that this assumes that 
the colour from the r-defined AUTO aperture would be the 
same as the colour from a r-defined Sersic aperture, how- 
ever, this is the closest estimation to a fixed Sersic aperture 
we can make at this time. 



7.3 Uncertainties within the photometry 

The gain value in SDSS data is constant within each stripe 
but varies between stripes. The SDSS mosaic creation 
process that is detailed in sections 13.21 and 13.41 combines 
images from a number of different stripes to generate 
the master mosaic. As the mosaics are transformed from 
different zeropoints, the relationship between electrons 
and pixel counts will be different for each image. This 
mosaic must suffer from variations in gain. The SExtractor 
utility can be set up to deal with this anomaly, by using 
the weightmaps generated by SWARP. However, this may 
introduce a level of surface brightness bias into the resulting 
catalogue that would be difficult to quantify. We calculate 
the SExtractor magnitude error via the first quartile value, 
taken from the distribution of gain parameters used to 
create the mosaic. The Gain used in the SDSS calculation is 
the average for the strip. The SExtractor error is calculated 
using Equation [S] where A as is the area of the aperture, a 
is the standard deviation in noise and F is the total flux 
within the aperture. By using the first quartile gain value, 



we may be slightly overestimating the component of 
the magnitude uncertainty calculation. However, given the 
amount of background noise in the mosaic, this component 
will constitute only a small fraction towards the error in the 
fainter galaxies, and in the brighter galaxies the uncertainty 
in magnitude due to the aperture definition will be much 
greater than the SExtractor magnitude error itself. The 
SExtractor magnitude error is calculated separately for 
each aperture type, and is available within the GAMA 
photometric catalogues. 



Am : 



F 
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We have attempted to quantify the uncertainty due to the 
aperture definition, in order to calculate its extent relative to 
the SExtractor magnitude error. We use the cleaned sam- 
ple defined in section |6~TI The dispersion in calculated mag- 
nitude between our different photometric methods for this 
sample are shown in Figures [9] 1101 111! [T2l and 1131 Figure 1211 
shows the relative scales of the uncertainty due to a galaxy's 
aperture definition (calculated from the standard deviation 
in AUTO/PETRO luminosities from the SDSS survey and 
our r/K / self-defined catalogues) and the error generated by 
SExtractor in the r band. The aperture definition uncer- 
tainty is generally much greater than that due to background 
variation and that SExtractor derives. Figure [22] shows 
how this standard deviation in a galaxy's r band magnitude 
changes with apparent magnitude. Whilst this uncertainty 
is larger than the SExtractor magnitude error, it is fun- 
damentally a more consistent judgement of the uncertainty 
in a given galaxy's brightness as it does not assume that 
any particular extended-source aperture definition is correct. 
Whilst the dispersion of the relationship increases with ap- 
parent magnitude (along with the number of galaxies), the 
modal standard deviation is approximately constant. Taking 
this to be a good estimate of the average uncertainty in the 
apparent magnitude of a galaxy within our sample, we have 
confidence in our published apparent magnitudes to within 
±0.03 mag in gri, ±0.06 mag in z, and ±0.20 mag in it. We 
calculate the same statistics in the NIR passbands (though 
without SDSS Petromag). We have confidence in our pub- 
lished apparent magnitudes within ±0.05 mag in YJHK; 
approximately two and a half times the size of the photo- 
metric rms error UKI DSS was designed to have (±0.02 mag, 
lLawrence et al.ll2007h . 



7.4 Number counts 

In order to construct a unbiased dataset, it is necessary for 
us to calculate the apparent magnitude where the GAMA 
sample ceases to be complete. 



7.4-1 Definition of GAMA galaxy sample used in this 
section 

The GAMA sample used in this section is defined as those 
SDSS objects that are within the area that has complete 
ugrizYJHK colour coverage, and have passed the star- 
galaxy separation criteria. Of the 908022 objects in the 
GAMA master catalogue, only 124622 fulfil this criteria. The 
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Figure 17. Sersic minus GAMA r-defined Auto magnitude against r-defined Auto magnitude, in all nine bands, for all objects in 
the GAMA sample that pass our star-galaxy separation criteria, and have credible ugrizY J H K r-defined Auto magnitudes. Contours 
increase geometrically in powers of 2, from 4 to 512. Bins are 0.1 mag (x axis) X 0.05 mag (y axis) in size. 



area of sky that has complete GAMA ugrizY J H K cover- 
age is 129.1232 ± 0.0008 sq deg; 89.7% of the entire GAMA 
region. All magnitudes in this section are r-defined AUTO 
magnitudes, unless otherwise defined. 



7.4-2 Determination of apparent magnitude limits 

Figure [23] shows how the sky density of GAMA galaxies in 
the nine passbands varies with apparent magnitude. The 
distributions peak in the 0.1 magnitude bins centred at 
w = 21.25, g = 20.55, r = 19.75, i = 19.25, z = 18.75, 
Y = 18.65, J = 18.45, H = 18.05 and K = 17.75 mag. Ta- 
bles[7][8]and[9]contain the number counts of GAMA galaxies 
in ugrizY J H K , this time using 0.25 magnitude bins. Pois- 
sonian uncertainties are also included. Both sets of data have 
been converted to deg -2 mag -1 units. 



The r band number count drop off, despite hitting the 
petromag_r= 19.8 mag GAMA main sample magnitude 
limit, is not absolute because the SDSS limit was extended 
to petromag_r= 20.5 mag in the GAMA 12 region so 
that useful filler objects could be selected, and because 
r&Aio/K/z band selected objects in G9 and G15 will also 
be included within the catalogue. Objects that are fainter 
than r mo[ j e ; = 20.5 mag (722 sources; 0.5% of the sample) 
will be due to differences in object extraction between SDSS 
and SExtractor, as mentioned in previous sections. 
The turnovers in Figure[23]will occur where the r = 19.8 mag 
limit is reached for galaxies with the median passband — r 
colour. We are within the domain where the number of 
galaxies within a magnitude bin increases linearly with in- 
creasing apparent magnitude, but a deviation from this re- 
lationship is visible in the figure approximately 3 magnitude 
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Figure 18. Sersic r - SDSS CMODEL r magnitude against Sersic index, SDSS r CMODEL magnitude, z, M r Sersie f° r objects in the 
GAMA sample that have passed our star-galaxy separation criteria, and have credible urK r-dcfincd Pctrosian magnitudes. Contours 
increase geometrically in powers of 2, from 4 to 512. 
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Figure 19. Sersic r - GAMA r-defined Petrosian r against Sersic index, r-defined r band Petrosian magnitude, z, M ri Sersic f° r au 
objects in the GAMA sample that have passed our star-galaxy separation criteria, and have credible urK r-defined Petrosian magnitudes. 
Contours increase geometrically in powers of 2, from 4 to 512. The brown function plott ed in the Sersic r - GAMA r-defined Petrosian 
r against Sersic index plot is taken from Figure. 2 (upper panel) of I Graham et all ll2005h . 
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stributions from GAMA r-defined Auto magnitudes, GAMA Sersic magnitudes and SDSS PETROMAGs (left to right), 
= u,g,i,z,Y ,J \H and K. Sources used are those within the subset region with credible Auto and Sersic magnitudes, 
DSS saturated object bit set. Contours increase geometrically in powers of 2, from 2 to 512 galaxies bin -1 . Bins are 
in size 
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Figure 21. The distribution of standard deviation in r band 
apparent magnitude against SExtractor's calculated magnitude 
error (using the first quartile gain from the gain distribution of 
the mosaic's input images) for our clean sample of galaxies, using 
SDSS, r-defined, X-defined and self-defined AUTO and PETRO 
magnitudes to calculate the standard deviation. Contours rise lin- 
early by 16 galaxies bin -1 , ranging from 8 to 120 galaxies bin -1 . 
Bins are 0.004 mag (x axis) X 0.001 mag (y axis) in size. 
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Figure 22. The distribution of standard deviation in r band ap- 
parent magnitude against apparent magnitude for our clean sam- 
ple of galaxies, using four different sets of magnitudes to calculate 
the standard deviation in each case. Contours rise linearly by 20 
galaxies bin -1 , ranging from 10 to 170 galaxies bin -1 . Bins are 
0.1 mag (x axis) X 0.01 mag (y axis) in size. 



Band 


Sources 


Redshifts 


% Redshifts 


Star-galaxy separa- 


124622 


82926 


66.5 


tion criteria only 








u<21.0 


46006 


39767 


86.4 


g^20.3 


67913 


58956 


86.8 


r^l9.8 


106032 


79672 


75.1 


K19.0 


74885 


66981 


89.4 


zsgl8.5 


59470 


55202 


92.8 


Ys;i8.4 


57739 


53339 


92.4 


J<18.2 


60213 


54264 


90.1 


HsS17.8 


55734 


51033 


91.6 


K<17.5 


46424 


43252 


93.2 



Table 6. The number of sources within the star-galaxy separa- 
tion and apparent magnitude limited GAMA samples that have 
a complete set of good ugrizY J H K r-defined magnitudes, the 
number of those sources that have redshifts from first and second 
year data and the percentage redshift completeness. Apparent 
magnitudes arc r-defined magnitudes, using the AB magnitude 
system. 



bins before the turnover occurs in all bands except r. This 
effect is due to colour incompleteness becoming a factor. Un- 
fortunately, despite our v&dio/K / z selection, there will be a 
population of objects that are bright in other passbands, but 
too faint in r to be included within our sample. Assuming the 
passband — r colour distribution is approximately Gaussian, 
this population will feature predominantly in the apparent 
magnitude bins near the turnover, causing the characteristic 
flattening we see. Accounting for this effect, we define the 
apparent magnitude sample limits of our sample to be a few 
bins brighter than this turnover, where the linear relation- 
ship still holds. Our apparent magnitude limits are set to 
u = 21.0, g = 20.3, r = 19.8, i = 19.0, z = 18.5, Y = 18.4, 
J = 18.2, H = 17.8 and K = 17.6 mag. 



7.4.3 GAMA apparent-magnitude limited catalogues 

Table [6] contains the sizes of the apparent magnitude lim- 
ited samples, and their current redshift completeness. Our 
magnitude limited optical samples contain approximately 
forty thousand less g alaxies than the equivalent samples in 
iBlanton et "ail (|2003l ). which covers the SDSS DR2 region, 
but extend two magnitude s deeper. When w e compare our 
number counts to 6dFGS l|jones et alj|2009l ); our samples 
are also smaller in area coverage, but si milar in s i ze an d 
much deeper in magnitude completeness. iHill et all lj2010( ). 
our previous attempt at defining a sample across the optical 
and NIR (combining MGC, SDSS and UKIDSS data to form 
a B band selected ugrizYJHK catalogue), was just one 
tenth of the size and was 0.2 — 1.8 mag shallower. Currently, 
a large fraction of our samples have not been spectroscop- 
ically sampled. After the completion of the 2008-2009 allo- 
cations of AAOmega spectroscopy, our apparent magnitude 
samples have ^ 75% completeness. It is anticipated that this 
statistic will have raised to > 95% after the completion of 
the 2010 allocation. 



7.5 Incorporating GALEX data 

We have also combined the GAMA sample with UV data. 
The GAMA master catalogue ha s been matched to GALEX 
photometry l|Wvder et al.ll2005l ). As GALEX observations 
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m r defined AUTO 

Figure 23. Number counts of GAMA galaxies (sources that have passed our star-galaxy separation criteria) with good ugrizY J H K 
colours, split into 0.1 magnitude bins and divided by the total area they cover. Error bars shown are for Poissonian number counts. 



are low resolution (typical imaging FWHM ~ 10 arcsec), the 
matching is complex compared to the simple UKIDSS/SDSS 
matching described within this paper as a number of sepa- 
rate SDSS objects may be matched to one larger GALEX ob- 
ject. The precise method of generating the GALEX matches 
is described in Robotham et al (in prep). In summary, all 
SDSS objects within the 90% Petrosian radius of a GALEX 
source are considered to be contributing flux to that source. 
The flux of the GALEX object is then apportioned to the 
SDSS objects, with the alloted fraction calculated via the 
distance between the SDSS and GALEX object. If no other 
nearby source is within 2.5 mag (in g) of the closest match, 
all flux is assigned to the closest match. GALEX has two 
distinct niters NUV and FUV. The generated magnitudes 
are stored within columns labelled MAG.AUTO-FUV and 
MAG-AUTO.NUV. 



7.6 SED fits using GAMA data 

The SEDs of 10 galaxies selected at random from the GAMA 
sample are shown in Figure [24] We show the GAMA-Galex 
UV luminosities, GAMA Total luminosities, Petrosian lu- 
minosities taken fr om the UKIDSS and SDS S surveys, and 
2 Bruzual-Charlot (|Bruzual fe Chariot! [2001 ) galaxy models 
with different ages and metallicities. The models are nor- 
malised via least squares best-fitting to the 9 GAMA data- 
points. For image clarity, we do not show the uncertainties 
on SDSS and UKIDSS datapoints. GAMA UV uncertain- 
ties are taken from SExtractor magnitude errors. GAMA 
optical and NIR uncertainties are calculated using the stan- 
dard deviation in the luminosity when different photometric 
methods are used (following the method described in Section 
I7.3[l . In some cases, the photometry provided by Survey- 
catalogues and the GAMA photometry are near identical, 
and match the galaxy models well (see 216774 and 137440 
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Figure 24. SEDs of 10 GAMA galaxies using GAMA matched GALEX photometry and r defined AUTO ugrizYJHK photometry 
(black triangles), and the comparable SDSS (blue circles) and UKIDSS (red circles) PETROMAG photometry. Uncertainties shown 
for GAMA ugrizY J H K points are calculated from the standard deviation in the photometry (as in Section I7.3|l . GAMA-GALEX 
uncertainties are SExtractor errors from the GALEX pipeline catalogues. Two Bruzual-Charlot 03 models are also plotted: the grey line 
is a 11 Gyr model using Z=0.05 and the purple line is a 0.25Gyr model using Z=0.02 (Z©). The models shown are the same as those in 
Figure [T6l 
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Table 7. Number counts for the ugr niters, using r-defined AUTO photometry, with Poissonian uncertainties. 



on Figure [24}. In other cases, where there is a discrepancy 
between our derived luminosities and the Survey-catalogue 
parameters, the GAMA photometry is a better fit to the 
models (e.g., 202588, 518102 on Figure EH)) . We therefore 
judge our r-defined AUTO colours to be a significant im- 
provement. 



THE IMPACT OF THE PHOTOMETRIC 
METHOD ON THE OBSERVED 
LUMINOSITY DISTRIBUTION 



8.1 Comparison between r band luminosity 
functions 



7.7 Released GAMA photometry 

The GAMA photometry described in this paper provides 
the GamaPhotom catalogue. This catalogue is filtered and 
combined with the other GAMA catalogues to produce the 
first GAMA data release, defined in Driver et al, 2010. 



In order to illustrate the effect that photometric methods 
have on statistical measurements of the galaxy population, 
we derive the r band luminosity function using 9 different 
photometric methods. We use all our photometric systems 
and the original SDSS photometry, and derive luminosity 
functions from the same population of galaxies. This should 
provide a consistent analysis for each method, removing all 
systematic effects except for that produced by the photo- 
metric method. 
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Table 8. Number counts for the izY filters, using r-defined AUTO photometry, with Poissonian uncertainties. 



8.2 Luminosity distribution and function 
measurement 

A number of techniques e xist for measu ring the galaxy 
luminosity distribution (see IWillmer] Il997fl , and functions 
to parameter i se it. We follow the methodology described 
in iHill et all l|20T0t ). This utilises the stepwise maximum- 
lik elihood method ( S WML ), originally described in detail 
in lEfstathiou et al.1 (|l98Sf ). and the sta ndard functiona l 
for m, the Sc h echte r luminosity function l|Schechterl Il97fj ). 
In IHill et all (|2010h . we developed a unique flux limit for 
each object based upon the spectroscopic limit and a colour 
limit. As we are working only in the r band, the colour limit 
is now used to calculate the change in the magnitude limit 
between the studied photometric method and SDSS petro- 
MAG. We now set the apparent magnitude threshold for each 
object using r ap erture, limit = 19.4 mag. An unfortunate side 
effect of the SWML method is that it requires normalisa- 
tion to calculate the luminosity density. We use the method 
of luminosity density scaling described in sections 3.1.2 and 



3.1.4 of lfiill et al.l (|2010T ). This involves calculating the num- 
ber density of galaxies within a 1 magnitude range contain- 
ing the M*mag galaxies, and using this to work out the 
required scaling multiplier. We also account for the cosmic 
variance within the GAMA regions. We calculate the source 
density of galaxies within a 5150 sq deg section of the SDSS 
survey (large enough for cosmic variance to be negligible) 
with dereddened -21.09 < M r - 5log 10 h < -20.09 (i.e. 
M* — 5logioh ± 0.5 mag, taking M* — 5logioh from the re- 
defined tauto photometry) and 0.023 < z < 0.1, and com- 
pare this with the source density calculated (using the same 
catalogue) from the GAMA regions of sky. We find that the 
GAMA regions are 95.2% as dense as the SDSS superpop- 
ulatior0, and therefore scale our 4>* parameters upwards by 

5 In essence, the GAMA survey is a post-stratified sampling of 
the SDSS, with the GAMA regions a stratum of the entire SDSS 
area. The SDSS source density is a universal parameter of our 
superpopulation, and can be used to improve the accuracy of the 
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488.99 ± 1.946 


17.875 


384.191 ±1.725 


17.875 


504.448 ± 1.977 
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Table 9. Number counts for the JHK niters, using r-dcfined AUTO photometry, with Poissonian uncertainties. 



a factor of Q * 52 . The area incompleteness of the K band- 
defined sample is accounted for by calculating the normali- 
sation volume with Area=133.5sq deg (the total coverage of 
the GAMA regions by K band UKIDSS data), rather than 
the Area=143.9sq deg used for the other samples. 

8.3 Sample selection 

We limit our sample using our star-galaxy separation cri- 
teria and an apparent magnitude limit of r ^ 19.4 mag 
(imposed on the dereddened magnitude system used to cal- 
culate the luminosity function). We use a brighter appar- 
ent magnitude cut than that defined in Section T7.4I because 
19.4 mag is the GAMA sample's target completeness limit 
over all three regions. Brighter than this limit our sam- 
ples are 91.3% spectroscopically complete (using tauto)- 
Our samples suffer greatly from spectroscopic incomplete- 
ness fainter than this magnitude limit. We impose a limit 

total luminosity density estimation we make from the GAMA 
dataset. 



based on the spectroscopic limit and a colour limit (i.e. 
19.4— (tsdss — r ap erture))- To remove the necessity of mod- 
elling the K or E corrections for each galaxy, we also impose 
a redshift limit of 0.0033 < z ^ 0.1. We adopt an evolution 
P = (where E(z) = 2.5,5 log 10 (l + z), setting (3 = de- 
notes no evolution in thi s redshift range) , and K(z) — 0.95z 
(following the r band in lHill et al.ll2010l ). We use the SDSS 
EXTINCTION_R parameter to deredden all our photomet- 
ric methods. We combine the data from the three GAMA 
regions, and treat them as one sample. Column 2 of Table 
1101 contains our sample sizes. 

8.4 The effects of surface brightness bias on the 
presented luminosity distributions 

Aperture selection can systematically bias the calculation 
of the luminosity distribution, particularly where a sample 
has a high surface brightness constraint (seelCross fc D river 
2002, particularly their Figure 5, and ICameron fc Driver 
[2007|). The SDSS photometric pipeline unfortunately is in- 
complete for Mr.sq > 23 mag arcsec -2 (see Section 3.4 of 
iBaldrv et al.l |2010| and references therein) . It follows that 
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0.02.x 
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1.99 
1.95 
1.98 
1.74 
1.90 
1.91 
1.97 
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0.11 
+0.13 
0.12 
+0.15 
0.13 
+0.14 
0.13 
+0.12 
0.11 
+0.15 
0.14 
+0.15 
0.14 
+0.14 
0.14 
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Table 10. The number of sources that pass our star-galaxy separation criteria, redshift limit and r ^ 19.4 mag limit, depending 
on which magnitude system is used to define the r band magnitude, with compariso n luminosity function parameters from SDSS 
llMontero-Dorta fc Pradall2009l , iBlanton et alj|2003h and SDSS+MGC llHill et aUhoioT) defined samples. All magnitudes use the AB 
magnitude syste m, and have been dereddened using the EXTINCTION_R SDSS parameter, j statistics are calculated using Mq r =4.71 
from Table 1 of iHill et al.l feOloT) . Note that the comparison stud y samples have much brighter magnitude limits; 17.77 mag in 
iMontero-Dorta fc Pradal , 17.79 mag in lBlanton et all and 18.76 mag in IHill et al.l 



any spectroscopic survey that bases itself upon SDSS pho- 
tometry, such as GA M A, will suffer from the same flaw. 
ICross fc Driver! |2002t ) have quantified the surface bright- 
ness dependency that the luminosity distribution inherently 
suffers from, and advise that a bivariate brightness distribu- 
tion (BBD) is the best way to quantify, and remove, SB bias. 
That is beyond the scope o f this paper, but this s hall be ex- 
plored in subsequent work. ICross fc Driverl (|2002T ) point out 
that a sample that is complete to fj,u m 24 mag arcsec -2 has 
very little uncertainty in its Schechter parameters due to SB 
selection effects, as the L* population th at define the fittin g 
are fully covered (see also section 4.1.2 of lDriver et al.| [2005). 
VST KIDS should provide such a catalogue. For now, how- 
ever, we accept that the SDSS input catalogue will not con- 
tain all faint, low surface brightness galaxies. The luminos- 
ity functions we present in this section are for samples that 
are surface brightness complete to fir.so < 23 mag arcsec -2 , 
and suffer from varying levels of completeness between 23 < 
/A-,50 < 26 mag arcsec -2 . As these luminosity functions are 
for a specifically low redshift sample, however, the effects of 
the surface brightness selection bias should be minimised. 

8.5 The effects of the aperture definition system 
on output Schechter parameters 

Figure [25] shows the luminosity distributions generated from 
different aperture systems, and illustrates how dependent 
the best-fitting luminosity function parameters are on the 
choice of aperture definition. The best fitting Schechter func- 
tion parameters (calculated via x 2 minimisation) are shown 
in Table ITOl 

The proximity of the r and self-defined mag luminosity dis- 
tributions signify that changing the SExtractor detection 
threshold (these catalogues utilise a detection threshold of 
1.7 a and 1 a respectively) has a limited effect on the proper- 
ties of a large sample (M* — blogioh ± 0.055 mag, a ± 0.014, 
4>* ± 0.0005 h 3 Mpc -3 ). There is an offset between the K 
and r-defined best fitting luminosity functions. This is not 
caused by the cosmic variance in the missing area of the 
K band sample; the best fitting Schechter function param- 
eters vary only slightly when this is accounted for. By us- 



ing the COVER-BITWISE flag, we can define a popula- 
tion of galaxies that are covered by K band imaging. The 
best fitting Schechter parameters for an rAuto sample within 
area covered by K band imaging (and normalised to the 
smaller volume) are M* — blogioh— -20.791 mag, a=-1.115, 
and 0*=O.O114 h 3 Mpc -3 - consistent with the area-complete 
LF within the uncertainty. It may be caused by a systematic 
alteration in the definition of the apertures used to calculate 
the flux of the galaxy population. 

The best-fitting elliptical Kron and Petrosian aperture lu- 
minosity functions are similarly distributed in the r, K and 
se£/-defined samples, indicating that the choice of light- 
distribution defined aperture does produce an offset that can 
be quantified. Whether the aperture is circular or elliptical 
is important. The SDSS petromag luminosity distribution 
should be similar to the r and self-defined elliptical PETRO 
distributions, but there is a noticeable M* — blogioh off- 
set (0.20 mag; inset of Figure l2"5fl . There is also a marked 
discrepancy between luminosity distributions calculated us- 
ing total magnitude apertures (Sersic and SDSS model- 
mag), and light-distribution defined apertures. The lumi- 
nosity distributions of the former are overdense for faint 
galaxies (M r — blogioh —16), and their best fitting power- 
law slopes are thus flatter. 

The Sersic luminosity distribution also measures higher den- 
sities of the brightest galaxies. This result is real, we have 
visually inspected the 139 galaxies that are distributed in 
the —22 ^ M r — blogioh ^ —21.5 magnitude bin and 
only 2 have suffered catastrophic failureqj. The remainder 
are generally well fit, though prominent spiral features do 
pose difficulties for the fitting algorithm. Of the 527 galax- 
ies with M T] sersic — blogioh < —21 mag within our red- 
shift limited, apparent magnitude cut sample, only 8 have 
m rt sDSS — Tn r , Sersic > 0. The marked discrepancy between 
the Sersic M* — blogioh parameter and that generated with 
the other samples (0.33 mag brighter) is indicative of a sce- 
nario where galaxies are moved out of the magnitude bins 
near M* — blogioh, and into the brighter bins. M* is not 

6 These profiles are viewable at 

http://star-www.st-and.ac.uk/~dth4/139cyc/ 
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Figure 25. Luminosity distributions, and the best fitting Schechter functions, calculated using different aperture definitions. Inset: 1 
sigma chi-squared best-fit contours in the M*-a plane. Errors on luminosity distribution points are Poissonian errors. 



an independent parameter, it is correlated with the other 
Schechter parameters, and accordingly the <j>* parameter has 
declined. The total luminosity density (j in Table \lO\i . whilst 
15% higher, is consistent with that generated by the SDSS 
model magnitude within uncertainties. As our Sersic mag- 
nitudes are not truncated, and the SDSS model magnitudes 
are truncated at 7R e for a de Vaucolouers profile / 3R e for 
an elliptical profile, the M* offset between these photometric 
systems is expected. We note that the residuals generated by 
this Sersic-fitting process expose the requirement of multi- 
component galaxy decomposition; many galaxies have some 
structure within their central bulges that the Sersic pipeline 
cannot model . A multi-componen t extension to this pipeline 
is detailed in iKelvin et all |20ld) , and that paper also ex- 
amines any discrepancies between the Sersic and multiple- 
component fitting results. No matter which aperture sys- 
tem we use, the luminosity distribution is overdense in the 



M r — 5logioh > — 16 magnitude bins when compared to the 
best fitting luminosity function. This indicates an upturn in 
the space density of galaxies at the dwarf-giant boundary, 
and the limitations of the single Schechter function fit. 
As noted in the introduction, the Schechter parameters gen- 
erated are for a sample that will suffer surface brightness in- 
completeness fainter than fj, = 23 mag arcsec -2 . In a future 
paper, we intend to account for this effect by undertaking 
a complete bivariate brightness analysis of the sample. The 
total luminosity densities we show here may therefore be 
systemati cally underdense due to surface brightness limita- 
tions (c.f. ICross fc Driver|l2002h . It is also apparent that the 
simple Schechter function parameterisation is no longer a 
good fit for the luminosity distribution of galaxies at fainter 
magnitudes; there is an obvious upturn in each sample that 
is not being modelled. As the Sersic photometric system is 
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the only system that accounts for missing light, it is the most 
effective way of calculating the luminosity distribution. 



9 SUMMARY 

The GAMA photometric pipeline has been designed to com- 
bine photometric data from a number of sources in a scien- 
tific and consistent manner. We have generated a series of 
large mosaics from imaging data taken from the SDSS and 
UKIDSS instruments using the SWARP utility, and under- 
taken optical and NIR photometry on an r and K band- 
defined sample of sources using SExtractor. We have also 
used a Galfit based pipeline (SIGMA) to generate Sersic 
ugrizYJHK magnitudes for all sources within the GAMA 
sample that pass our star-galaxy separation criteria. 
We have created a set of r-defined source catalogues in 
the ugrizY J H K passbands, calculated the apparent mag- 
nitude limits at which these samples are complete, and es- 
timated their current redshift completeness. Whilst these 
samples do not contain as many sources as those used 
to ca l culate the SDSS luminosity fu nctions (|Blanton et al.l 
120031 . iMontero-Dorta fc Praddl2009h . they are complete to 
a greater depth. Our NIR apparent magnitude limited cat- 
alogues are of comparable size to those produced by the 
6dFGS ones et all [2009), with the decrease in coverage 
area matched by increasing depth. Our source catalogues 
have the advantage that they can be used to accurately com- 
pare the optical and NIR luminosity of any source within 
our sample. As the aperture definition is constant we can 
calculate exact colours. We have attempted to quantify the 
level of uncertainty in source flux due to the definition of the 
aperture. We find that the uncertainty is generally ~ 400% 
of the SExtractor output uncertainty (produced by back- 
ground noise variations and -ft). 

We have attempted to quantify the percentage of mis- 
matches between our catalogues and the SDSS catalogues. 
We have taken a subset region of our images and tested 
their results against the SDSS detection software. We have 
detailed reasons for mismatches between these catalogues. 
They are generally an issue of different deblending decisions, 
artifacts, background noise or proximity to saturated stars. 
We find that there do seem to be some low surface bright- 
ness, compact objects that are missed by the SDSS detection 
software that are found by SExtractor, and some false SDSS 
detections. These faulty sources would be removed from the 
spectroscopic observing list by the extensive visual classifi- 
cation process undertaken by the GAMA team. 
As our Sersic magnitudes are generated from the images that 
the profile is calculated on, the colours produced from the 
Sersic magnitudes are liable to be inaccurate and suffer from 
bias due to the aperture definition. We have defined a set of 
'Fixed aperture' Sersic magnitudes using the Sersic r band 
magnitude and colours from the r-defined AUTO catalogue. 
We recommend the use of these magnitudes, called X to tai- 
These magnitudes will provide both an accurate estimate of 
the luminosity of each galaxy, and optimal colours. 
The GAMA photometric pipeline catalogue does not have 
the breadth of focus that the existing products from the 
SDSS and UKIDSS surveys have. As both the instrument 
gain (in SDSS only) and seeing change from frame to frame, 
our mosaics have variations within them that would cause 



errors in the calculation of stellar photometry. We have not 
attempted to deal with saturated objects, which causes un- 
certainty at the bright end of our catalogue, or masking of 
artifacts. These problems primarily affect stars, or are re- 
moved following the extensive visual classification program 
we have undertaken. For the task our catalogues were devel- 
oped for, the accurate calculation of extended source aper- 
tures across instruments for a specific sample of extended 
sources, they are more adept than the existing tools. 
Finally, we have generated r band luminosity distributions 
and best-fitting luminosity functions from our population 
of galaxies using 9 different aperture definitions: SDSS cir- 
cular Petrosian and Model magnitudes, r-defined elliptical 
Petrosian and Auto magnitudes, A"-defined elliptical Pet- 
rosian and Auto magnitudes, sel /-defined elliptical Pet- 
rosian and Auto magnitudes, and an elliptical Sersic total 
magnitude. We find that the r and seZ/-defined elliptical, 
light-distribution defined apertures produce similar results, 
indicating that the choice of detection threshold is unimpor- 
tant. We find that there is a similar AUTO-PETRO offset 
in the r, self and A"-defined samples. We find that the use 
of circular apertures does have an effect on the best-fitting 
Schechter fit, with the SDSS petromag having a fainter 
M* parameter than the SExtractor samples. We also find 
that the use of total magnitude systems affects the slope 
of the luminosity function, with both the Sersic and SDSS 
modelmag luminosity functions having a steeper a parame- 
ter. When we calculate the total luminosity density for each 
sample we find that using the Sersic magnitude system gives 
us a higher value, approximately 15% higher than samples 
that use other aperture definition systems. Following visual 
classification of a subsection of our sample, it is clear that 
this is not due to errors within our Sersic magnitude calcu- 
lation. We also note that the Schechter luminosity function 
does not provide a good fit at the faint end of the luminosity 
distribution, and a clear upturn at the dwarf-giant boundary 
is seen. 
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The UKIDSS project is defined in lLawrence et al.ll2007l . 
UKIDSS uses the UKIRT Wide Field Camera (WFCAM; 
Casali et al. 20 071) and a photometric system described in 
Hewett et al. | l200d The pipeline processing and science 
arc hive are described in Irwin et al (2010, in preparation) 
and lHamblv et ail |200S| . 

STILTS is a set of command-line based catalogue 
matching tools, based upon the STIL library. It is currently 
being supported by its author (Mark Taylor), and can be 



downloaded from http://www.star.bris.ac.uk/~mbt/stilts/ 
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Figure Al. The effects of convolution and the change in passband of observations of SDSS object 588848900968480848. 
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APPENDIX A: VARIATION BETWEEN 
PASSBANDS 

Figure lAll shows the 18 200x200 pixel images of the piece of 
sky containing SDSS object 588848900968480848; 9 cutouts 
from the standard image mosaics, and 9 from the convolved 
image mosaics. What is easily noticeable is that not only 
does the ability to see features of the object change dra- 
matically between the u (top left) and K (bottom right) 
wavebands (spiral arms are visible in the optical, but in the 
K band there only seems to be a bar and a bulge compo- 
nent), but that objects around it appear and disappear (a 
small blip to the SE in the r band that may or may not be 
part of the object itself, at least 5 faint objects in the E of 
the frame in the NIR) . The size of the object seems to halve 
from the g band to the J band, though this may be an effect 
of the image quality (the SDSS g band should have a much 
smoother background than the UKIDSS J). The apparent 
magnitude of the object itself changes by 2.8 magnitudes 
from u to its peak in H (SExtractor calculates AB mag- 
nitudes using an r band-defined AUTO aperture of 16.67, 
15.44, 14.85, 14.49, 14.32, 14.23, 14.04, 13.83, 14. 11 mag in 
ugrizYJHK). This is probably due to the decrease in dust 
opacity from the UV to the NIR. 

The convolved images also show greater variation between 
the object and the background (these images all use a lin- 
ear scale between the 99.5% quantile pixel and 0, the back- 
ground). For instance, the extended spiral arm to the left of 
the bulge in the u band becomes slightly more apparent in 
the convolved u band image. The size of the object in the 
convolved images generally looks larger than the standard 
images, though this again is probably due to the smoothing 
of the background making flux overdensities more apparent 
in the convolved images. 



